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SECTION I 
INTRODUCTION 



Many problems encountered in scientific or engineering work involve mathematical 
equations or sets of equations whose solution in most cases is difficult or practi- 
cally impossible to obtain by the classical approach to equation solution. The TR-4-8 
Analog Computer provides the technical worker with a general purpose computer which 
permits the rapid solution of .linear or non-linear equations. 

Although the analog machine is correctly termed a computer, it does not perform its 
computations by serial calculations as does the desk calculator or digital computer. 
Instead it performs the required mathematical operations in a parallel manner on 
continuous variables. In the TR-4-&, as in most modern analog computers, the contin- 
uous variables are direct current voltage s . The electronic analog computer makes it 
possible to build an electrical model of a physical system where the voltages on the 
computer represent the dependent variables of the physical system. Except for a 
constant of proportionality, or scale factor, each voltage will behave with time in 
a manner similar to the physical system variable. Thus, if the vertical position of 
the center of gravity of an automobile oscillates with time during a disturbance, 
then the voltage representing the height of the center of gravity above the road sur- 
face will also oscillate j if the temperature of the coolant at the exhaust port of a 
condenser rises exponentially to a steady value, then so will the voltage representing 
it on the computer. 

It can be said that the actual system and the electrical model are analogous in that 
the variables which demonstrate their characteristics are described by relations 
which are mathematically equivalent. The actual system has thus been simulated be- 
cause of the similarity of operation of the electrical model and the physical system. 
This capability of the analog computer is of great value in performing scientific re- 
search or engineering design calculations because it permits an insight into the re- 
lationship between the mathematical equations and the response of the physical sys- 
tem. Once the electrical model is completed, well-controlled experiments can be per- 
formed quickly, inexpensively, and with great flexibility to predict the behavior 
of the primary physical system. 

Although the analog computer utilizes electronic components in its operation, it is 
not essential that the user have an extensive knowledge of electrical circuits. The 
TR-4-B is basically a set of mathematical building blocks, each able to perform spec- 
ific mathematical operations on direct voltages and capable of being easily inter- 
connected. By appropriately interconnecting these building blocks, an electrical 
model is produced in which the voltages at the outputs of the blocks obey the rela- 
tions given in the mathematical description of a physical problem. 
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Since our interest is frequently in the dynamic behavior of physical systems, the 
isatheiaatical equations are usually differential equations having time as the inde- 
pendent variable,, In order to solve such equations, the standard components of the 
computer must perform the following operations? inversion, algebraic summation, in- 
tegration with respect to time, multiplication and division, and function generation. 

The sequence of steps for constructing a dynamic model on an analog computer requires 
first a mathematical description of the physical system, usually in equation form. 
From this description the operator derives the information necessary to set up a 
computer program for interconnecting the computing components and determines the 
required initial conditions and forcing functions . The computing components are 
interconnected with wires called patch cords. The input and output terminations of 
the computing components are brought out to a patch bay which is fitted with a re- 
movable patch panel into which the patch cords are inserted. The problem patching 
may, therefore, be done in advance, away from the computer. The problem is placed 
on the computer by inserting the patch panel and adjusting the problem parameters 
to the value of the first case to be investigated. Selected voltages are applied 
to various components in the form of inputs or initial conditions. These voltages 
are derived from a precise reference voltage. 

Once the computing elements have been patched, adjusted, and energized, the computer 
is switched into the operate mode. The voltages on the computer change with time 
in accordance with the equations that govern the physical system variables. The be- 
havior of the computer model is viewed through an output device such as an X-I plot- 
ter, oscilloscope, strip-chart recorder, or digital voltmeter. 

The intention of this manual is to provide the scientist or engineer, using an ana- 
log computer for the first time, with an introduction to its functions, programming, 
and operation so that he is able to achieve usable experimental results. The 
reader ought to gain from the early sections sufficient knowledge of the computer 
operation to enable him to interconnect computing components and to operate the com- 
puter without difficulty. Later sections provide an understanding of simple pro- 
gramming procedures and computing techniques so that he is able to construct a com- 
puter program for any straightforward investigation. Needless to say, the full story 
cannot be told in a short manual and in many cases the literature cited in the bib- 
liography provides a more comprehensive treatment. However, the ideas and facts pre- 
sented here should set the reader well along the best path towards dynamic problem 
investigation by effective use of the modern general purpose analog computer. 
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SECTION II 
THE TR-48 COMPUTER 



1. GENERAL DESCRIPTION OF THE TR-48 

The PACE® TR-48 (Figure 2.1-1) is a fully transistorised, general purpose analog 
computer. Consisting entirely of solid-state circuit elements, the TR-4.8 is com- 
pact in sise and is suitable for desk top mounting. The computer is able to operate 
stably and accurately in normal office ambient conditions,; there is no need for 
large primary power systems or special cooling duct installations since the power 
requirements are small. Reliable, with simplicity in functional design, the TR-48 
is easy to use and can be a powerful aid to the individual engineer in the rapid 
solution of scientific and engineering problems. 

The TR-4.8 utilizes a building block concept , in which individual computing compon- 
ents may be easily interconnected to solve the required equations by forming elec- 
trical models analogous to the system under study. Each building block, either 
individually or in combination with others, is capable of performing one or more 
of the following operations on variable DC voltages s multiplication by a constant, 
algebraic summation, integration with respect to time, multiplication of two vari- 
ables, and generation of known functions of a variable. Each component has input 
and output terminations which are readily accessible at the computer FTe-Patch 
Panel for interconnection by bottle plugs (jumpers) and patch cords. 

The Pre-Patch Panel is arranged in a series of twelve similar modules , with each 
module terminating a complete set of computing components. (Figure 2.1-2.) The 
modular design tends to eliminate patch-panel clutter caused by long across-the- 
panel patching. In addition, problem patching, checking., and trouble-shooting are 
more readily accomplished, and patching errors are less likely. The five patching 
areas within each module are color coded in vivid contrasting colors with large, 
clear lettering to help in faster and surer patching. The Pre-Patch Panel is re- 
movable to permit problem storage and also problem patching without committing 
operating time. 

To the left of the Pre-Patch Panel is the monitoring and control panel. This area 
contains the control switches, and metering circuits which permit turning the com- 
puter on and off, engaging and dis-eng aging the removable Pre-Patch Panel, and con- 
trolling the computer mode,, Monitoring facilities consist of a digital voltmeter 
and a multi-range voltmeter. The di gital voltmeter permi ts setting__pr oblem para - 
meter s to two-decima j L _ac^uracy and observing or measuring a lX""c omput ati oral and, 
operating potentials. 

To the right of the Pre-Patch Panel is the coeffieiesrk-setting potentiometer and 
function switch panel. This panel provides space for a maximum of sixty potentio- 
meters and five function switches. Fou r of every five potentiometers have o ne 
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end grounded while the fifth potentiometer is iiagrounded. The potentiometers and 
function switch connections are terminated in the appropriate areas as indicated on 
the Pre-Patch Panel. 

2. OPERATING CONSIDERATIONS 

The TR-4B is shipped complete with all components in place . The unit is completely 
calibrated and adjusted at the time of manufacture ♦ After performing the simple 
installation checkout procedure outlined in the maintenance manual, and connecting 
the unit to a suitable power source, the computer is ready for operation. 

The use of low level reference supplies (+10 and -10 volts) with their associated 
current-limiting circuits eliminates shock hazards when patching with the Pre- 
Patch Panel in place in the computer. The current-limiting circuits also protect 
the reference supplies from damage due to skorting to ground, or to each other ♦ 
Thus an errant patching connection (shorting the plus reference to the minus ref- 
erence for example) will not adversely effect the supplies (output current drops 
to aero) nor will the reference supply fuses blow. In addition, the metal portions 
of the Pre-Patch Panel (except the handles) have a scratch-resistant, non-conduc- 
tive paint coating which, in conjunction with the plastic patching blocks, gracti- 
c &£ly eliminates shorting-out hanging patch cords. ^ <° ° <} 

a.. Preliminary Operating Considerations ^ """" 

The following steps are recommended prior to operating the TR-48 to prevent possi- 
ble false trouble indications. 

(1) Ascertain that gJJ operational amplifiers have four- connect or bottle- 
plugs properly placed and seated as shown in Figure 2.2-1 . This provides the ampli- 
fiers with feedback and prevents them from overloading during the problem solution. 
The procedure for inserting and removing the Pre-Patch Panel is described in Sub- 
Paragraph Jio 

(2) Patch the Digital Voltmeter (INK) and multi-range voltmeter (VM) to 
the selector readout system. (See Figure 2.3-2 and Paragraph 3j&« of this section 
for a description of the selector system operation.) 

(3) Apply power to the computer and depress the PS (potset) button. Ini- 
tally the overload lamps of the operational amplifiers will light j after a few 
seconds all of the lamps should extinguish. 

(4) Check the various supply voltages of the TR-i48. All power supply 
outputs are connected directly to the voltmeter FUNCTION switch (through appro- 
priate scaling resistors); thus, the check may be accomplished simply and rapidly. 
(See Paragraph 3s.. of this section.) 

(5) Check the plus and minus reference supplied for readout on the DVM by 
selecting M9 and A50. (See Paragraph 3a*) 
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FIGURE 2.1-1 TYPICAL TR-48, FRONT VIEW 
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FIGURE 2.1-2 PRE-PATCH PANEL MODULAR LAYOUT 
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FIGURE 2.2-1 AMPLIFIER WITH FOUR-CONNECTOR ' BOTTLE-PLUG PROVIDING FEEDBACK 



CAUTION 



If the VM FUNCTION switch is in the PATCH 
position the selected voltages are also 
applied to the VM. Be sure the RANGE 
switch is set to a high enough range selec- 
tion to handle these inputs * 

(6) Allow a few minutes warm-up time; this assures that the computing 
components , (including the DVM) are up to normal operating temperature. Ground the 
DVM input termination (designated DVM on the 12*763 READOUT PANEL) and adjust the 
DVM zero control on the back of the Control Panel door ( Figure 2.2-2 ) for a 
00,00 reading o Should the polarity relay begin to chatter, turn the zero control 
slightly clockwise until the chatter stops and the indicators retain the 00.00 
display. The DVM zero adjustment should be checked daily; initially, due to aging 
of the precision components, this adjustment may be required more frequently. 

(7) In the pot-set mode of the computer (PS button on the Control Panel 
depressed), closed relay contacts provide a feedback circuit for the operational 
amplifier. (See Paragraph U of this section for a more detailed description.) 
This feature permits the removal of the Pre-Patch Panel to balance the operational 
amplifier. However, when the computer is switched from pot-set to any ether mode, 
the relay contacts open and the circuit as patched on the Pre-Patch Panel provides 
the feedback loop. Momentary amplifier overload may result during the relay oper- 
ating time; thus to eliminate possible error in the computer solution, the operator 
should always switch the computer to re-set (depress the RS button) before switch- 
ing to the operate mode. This permits the pot set relays to open and the ampli- 
fier summing junctions to settle before starting the problem solution. There is no 
actual aa^^ag period required^ that is, the operator may depress the RS button 
and then immediately depress the button for the desired mode. This sequence of 
operation will prevent the possible momentary overloads from effecting the problem 
solution. 

Pre-Patch Panel Insertion and Removal 

To insert the Pre-Patch Panel, set the lip on the lower edge of the panel in the 
guide groove (Figure 2.2-3) <» Push the top of the panel in so the panel is verti- 
cal. (A micro-switch will not permit the Pre-Patch Panel motor to function un- 
less the panel is properly seated.) Applying a slight hand pressure to the cen- 
ter of the panel (to maintain the panel in the vertical position) depress and 
hold the ENGAGE button on the Control Panel. (The computer must be in the pot- 
set mode.) This button properly seats and firmly holds the Pre-Patch Panel in 
position. Note that the ENGAGE and DISENGAGE positions of the switch are spring 
loaded; therefore, the switch must be depressed and ,foe]4 until the limit switches 
stop the motor operation. This feature eliminates the possibility of accident- 
ally engaging or disengaging the Pre-Patch Panel. 

The removal of the Pre-Patch Panel is essentially the opposite of inserting the 
panel. When disengaging the panel (depress and the DISENGAGE button) a light 
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hand pressure should again be placed on the panel to eliminate the possibility of 
the panel falling for -ward when disengaged* 

£<> Amplifier Balance 

The d-c operational amplifiers are chopper stabilized to prevent drift and resultant 
errors in the computer results * Grift in an amplifier results in an output voltage 
(or s££&Sk) vith a zero input <> To eliminate offset, the amplifiers of the TR-48 
are balanced , i.e., vith a zero input, a bias current is applied to the amplifier 
summing junction equal and opposite to any current due to drift thus placing the 
summing junction at virtual ground. Once balanced, drift in the amplifiers is 
eliminated automatically by the stabilizer circuit. 

The 6>c amplifiers of the TR-48 are extremely stable and normally do not requir e 
b alancing for periods up to several months . To assure accuracy and confidence 
in the computer results, it may be desirable to check the amplifier balance daily | 
this check can be made rapidly and simply since the selector system and voltmeter 
are used. The following is a step by step procedure for checking amplifier balance. 

(1) Place the voltmeter FUNCTION switch in the BAL position and depress 
the PS push-button of the MODE switch. 

(2) Using the selector system, select each amplifier ADO through A49 (A4& 
and A49 are the plus reference and minus reference amplifiers respectively). The 
voltmeter should register a zero deflection for each amplifier. 

(3) Should an amplifier cause a deflection to either side of the center 
zero on the voltmeter, adjust the corresponding balance control. The balance con- 
trols for amplifiers A00 through A47 (the operational amplifiers) are located di- 
rectly behind the Pre-Patch Panel (Figure 2.2-4a) . The balance controls for A4& 
and A49 (reference amplifiers) are located on the Reference Regulator 43.104 be- 
hind the Potentiometer Panel (Figure 2.2-4b). Adjust these controls for a zero 
reading on the meter. 

£0 Changing Computational Components 

In the solution of some problems it may be necessary to add a specific type of 
computational component to the existing complement. Since many of the module po- 
sitions are designed to handle more than one type of computing component, a com- 
ponent not required in the problem, investigation may be removed and another unit 
placed in that cradle. Figure 2.2-5 illustrates the various positions of the com- 
puting components in the TR-4S module area. This diagram illustrates which type 
of computing component is compatible with each cradle or module position. The 
procedure for replacing a computer component and changing the Pre -Patch Panel 
patching block is described in the following Sub-Paragraphs (1) and (2). 
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FIGURE 2.2-2. DVM ZERO ADJUSTMENT LOCATION 
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FIGURE 2.2-3 PRE- PATCH PANEL INSERTION 
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FIG'JRF. 2.2-4. AMPLIFIER BALANCE CONTROL LOCATION 
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MOTE 



Failure to change the Pre-Patch Panel 
patching block may prevent proper use 
of a computer component due to the 
arrangement of jumpers on the rear of 
the patching block o 

(1) Computing Module Replacement 

(aj Remove the Pre-Patch Panel to expose the component modules . Re- 
move the two phillips-head retaining screws from the top and bot- 
tom of the module to fee removed (Figure 2.2 -6). 

(Jk) Insert the special module removal handle in the holes provided 

in the central area of the module (Figure 2*2-6). Pull the module 
forward removing it from the TR~48 0 

(Si) Place the new component in place ? be sure the guide pins are pro- 
perly seated in the guide~pin holes before mating the connectors 
at the rear of the module , 

(j|) Check that the module is properly installed (connector firmly 
mated, etc.) and replace the tw retaining screws <> 

(2) Patching Block Replacement 

(aj The Patching blocks of the computing components are held securely 
in place by the retaining strips on the front of the Pre-Patch 
Panel (Figure 2.2-7a). 

(&) The retaining strip above or below the patching block may be 
removed to change blocks • The retaining strip is released by 
removing the four screws directly behind the strip on the rear 
of the Pre-Patch Panel (Figure 2o2-7b) . 

(&,) Once the retaining strip is free, remove the original patching 
block and replace it with the new block (Figure 2.2-7) . Secure 
the retaining strip with the four screws. The Pre-Patch Panel 
is now ready for problem patching. 

NOTE 

If blocks on two adjacent horizontal rows 
are to be replaced it is only necessary to 
remove the retaining strip between the two 
rows. Patching blocks may then be removed 
from the row above and below the strip. 
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3. MONITORING AND CONTROL 



The control panel of the TR-48 is designed to allow simple control and monitoring 
of the computer components (Figure 2.3-1) • The following sub-paragraph describes 
the function and operation of the various switches and controls mounted on this 
panel. 

ji. Signal Selector 

The signal selector consists of three vertical rows of pushbuttons s The first row 
contains two buttons designated A and P, and the second (tens) and third (units) 
rows contain buttons designated 0 through 9. Depressing the A button permits the 
operator to select the outputs of the 48 operational amplifiers (A00 through M7), 
the two reference supply amplifiers A48 and M9 (plus and minus reference respec- 
tively), and 15 of the IN trunk lines (selected as A$0 through A64). Depressing the 
P button permits the operator to select the individual wiper outputs of the sixty 
potentiometers (POO through P59). 

Depressing the tens button (assuming the A button is depressed) sets the selector 
system for amplifiers in a given tens group, i.e., if the 2 button in the tens row 
is depressed, amplifiers A20 through A29 are set up for selection. The units but- 
ton determines which of these ten amplifiers is actually selected. The selector 
system button numbering corresponds with the amplifier and potentiometer designa- 
tions as marked on the Pre-Patch Panel and Potentiometer Panel respectively. 

Once the letter and tens selector buttons are depressed, say P1 (addressing P10 
through P19), the individual potentiometers of this group may be selected by 
merely changing the units designation button, i.e., if P10 is selected, P11 may 
be selected by depressing the units 1 button only. 

The selector system output is connected to the three terminations marked SEL in 
the upper white portion of the 12.763 Readout Panel (Figure 2.3-2a) . 

In order to read out a selected signal on the DVM, the upper bottle plug shown in 
Figure 2.3-2b must be in place. The multi-range voltmeter (VM) may also be con- 
nected to the selector line by installing the lower bottle plug shown in the illus- 
tration (the VM FUNCTION switch must be in the PATCH position) . It should be noted, 
however, that the VM circuit will load the output of the selected component with a 
relatively low impedance and should not be used if the monitored circuits cannot 
tolerate this load. (Set the FUNCTION switch to a position other than PATCH to 
disconnect the VM from the SEL terminations.) The VM should not be used when set- 
ting attenuators . To prevent this possibility from accident ly occurring, a relay 
circuit disconnects the VM from the selector system when the computer is ^placed in 
the pot— set mode. 

Digital Voltmeter 

The Transistorized Digital Voltmeter 26.183 is terminated in the 12.763 Readout 
Panel area (Figure 2.3~2a) and is designated DVM. As previously mentioned the 
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FIGURE 2.2-6 REMOVAL OF COMPUTING MODULE 
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FIGURE 2.2-7 PATCHING BLOCK REPLACEMENT 
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FIGURE 2.3-1 TR-48 CONTROL PANEL 
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READOUT PANEL 
12.763 



READOUT PANEL 

12.763 



FIGURE 2.3-2 READOUT PANEL 12.763 AREA 
AND DVM/VM TO SEL PATCHING 



DVM may be bottled to the SEL termination to monitor the selector system voltages, 
or -with a patch cord the DVM may be used to monitor signal levels at practically 
any termination in the Pre-Patch Panel , The WM has a 10 megohm (minimum) full- 
time input impedance. 

for -the. JVM 

The only operational setting required. is periodic adjustment of the zero control* 
(See Paragraph 2&,, Operating Considerations o ) 

Multi-Range Voltmeter 

The Multi-Range Voltmeter is permanently wired into various circuits of the TR-48 
to facilitate rapid readout of certain voltages by selection with the meter FUNC- 
TION switch. The voltmeter also has a RANGE switch, thus permitting close to full 
scale readouts for maximum accuracy. The ranges are 0,1, 0.3, 1,3, 10, and 30 
volts, (in addition to an off position). The RANGE switch functions only when the 
function switch is in the PATCH position. The PATCH position connects the volt- 
meter to the Pre-Patch Panel VM termination (Figure 2,3-2) permitting this point 
to be bottle-plugged to the SEL output or, as in the case of the WM, monitoring 
voltages at most Pre-Patch Panel terminations via a patch cord. 

Following is a list of the voltmeter FUNCTION switch positions with a brief des- 
cription of the function of each position. 

BAL 



PATCH 



RELAY 



-15, -3, 30, 15, and 2 



Computer MODE 

The operating mode of the computer is controlled by the six-pushbutton selector 
just to the right of the POWER button (Figure 2.3-1 )» Following is a list of the 
pushbutton and brief description of their functions, 

MODE BUTTON DESCRIPTION 

OP (Operate) When this pushbutton is depressed, all integra- 

tors are simultaneously released to respond to 



DESCRIPTION 

Connects the stabilizer output of the amplifier 
addressed by the selector system to the meter 
to facilitate checking and/or adjusting the 
amplifier balance. 

Connects the meter, via the RANGE switch, to the 
VM Pre-Patch Panel termination. 

Connects meter to the relay power supply (-20 
volts) . 

Connects voltmeter to output of corresponding 
power supply. z^&C ci^on^cJ^^J^ <zjlJ.&&Z*l 
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OP (Operate) (Continued) 
HD (Hold) 



RS (Reset) 



PS (Pot Set) 



SL (Slave) 



RO (Rep-Op) 



S.c Overload Indicators 



input signal voltages* The integrator outputs change 
in potential as dictated by the inputs; a time vary- 
ing behavior is produced* This generates the voltage 
solution of the programmed problem* 

Depressing the RD pushbutton permits the problem solu- 
tion to be stopped and all voltages held at the poten- 
tial attained up to the instant of depressing the but- 
ton* The problem may be continued from this point by 
depressing the OP button or re-set to the starting 
point by depressing the RS button* 

In the RESET mode all circuits except the integrators 
function normally* The integrator outputs are held 
at their respective initial conditions (IC) as dic- 
tated by the IC input voltage* (The integrator out- 
put is zero if no IC voltage is applied*) 

Amplifier input resistor summing junction groundedj 
permits setting potentiometers under actual load* Al- 
so provides amplifiers with relay-contact feedback so 
Patch Panel may be removed to balance amplifiers* 

When a TR-48 Computer is to be slaved to another TR-4.8 
(master) this button is depressed. The slaved computer 
then responds to the selected modes of the master com- 
puter pushbuttons. 

The RO button switches the computer into the repetitive 
operation mode* The computer switches automatically 
between operate and reset at a predetermined rate, 
(See Paragraph ^ of this section*) 



The overload indicators (Figure 2*3-1) provide a visual alarm when an overload oc- 
curs in any of the operational amplifiers, i*e*, when the summing junction is not 
at virtual ground. An overload may be due to improper scaling, improper patching 
or loading. 

When the computer is initially turned on, all the indicator lamps may light,; however, 
in a few seconds, as the amplifiers settle, all the lamps should go out. Should a 
lamp remain lit it could be caused by a patching error such as the failure to pro- 
vide un-used amplifiers with feedback (via the four connector bottle plugs). Pro- 
longed overloads will not damage an amplifier. 
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£. Trunks 



The trunks (terminating at the Trunk 12.762 area) provide point-* o-point connections 
to the connectors at the left-rear of the computer (Figure 2.3-3) . These connectors 
may be used as outputs to accessory equipment, or the trunk terminations may be 
cabled to a second TR-48 as signal carrying lines for the interconnection of the 
problems patched on the separate Pre-Patch Panels of slaved computers. 

g,. Readout Devices 

The problem solution obtained from the TR-48 may be permanently recorded or tempor- 
arily displayed on various types of readout devices. A few of the more common read- 
out devices are described in this paragraph. 

Either the WM or voltmeter may be used to read out the problem solution. These 
devices, however, do not record the solution and thus a permanent record is not 
made. When in REP-OP the computer solution is, continuously displayed on an oscillo- 
scope; in this case a permanent record may be obtained by photographing the scope 
trace. X-Y plotters (such as the EAI 1110 VAR IPLOTTER® ) or strip-chart recorders 
may be used for a permanent record of TR-48 problem solutions. These units, how- 
ever, do not have the frequency response necessary to accurately record the solution 
when the computer is placed in the high-speed REP-OP mode of operation. The Readout 
Panel scope and plotter terminations are wired to connector plugs at the rear of the 
TR-4#. See Figure 2.3-3. 

4. ATTENUATORS 

One of the simplest and most useful operations performed on an analog computer is 
the multiplication of a variable voltage by a positive constant less than unity; 
i.e., attenuation of a signal- The TR-48 has a basic complement of 10 attenuators 
(or potentiometers) and may be expanded to a full complement of 60 potentiometers. 

Each Attenuator Group 4.2.283 provides five potentiometers for setting problem co- 
efficients, initial conditions, and problem inputs. The potentiometers are mounted 
in up to 12 horizontal rows of five potentiometers per row. Each row is terminated 
at an individual patching area of the Fre-Patch Panel; space is provided for one 
potentiometer patching area per module of the Pre-Patch Panel. Four of the five 
potentiometers have one end grounded while the fifth has both ends ungrounded. (See 
Figure 2.4-1 «) 

The standard potentiometers in the TR-48 are 10 turn, wirewound, 5000 ohms , individ- 
ually fused units with calibrated dials and a locking mechanism. ~ 

The p otentiometer may be used i nconj unction with reference to obtain a fix ed 
laccurate volta ge less than ref erence, or_j^o_jauIt^ by any con - 

st ant less than janity. Figure 2.4-2 is a schematic of~~a potentiometer with +10 
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volts applied to the high end*; the output at the viper is k times +10, -where fc 
is: 

k - ^1 (EQ. 2.4-1) 

The potentiometer shown in Figure 2,4-2 is unloaded, and the mechanical ratio of 
Ei:% equals the electrical ratio e 0 :ei n ; thus, the potentiometer may be set to 
the exact ratio by means of the calibrated dial attached to the wiper shaft. How- 
ever, the two ratios will not be equal when the potentiometer is loaded as is the 
case when it is used as a computer problem element. Normally, the pot is loaded 
by^ithej^^JOO^QOO ^QQK) or 10,000 (1Q K)_ohms resistor since _ a~potent iometer ^n- 
erally_j ^ds^^ampljfier~ and these values are the most commo n ~amplifier~ ~inpu t~re ; ^ 
sistors_, Figure~2.4-3 illustrates tne errect on ^he~e±^ie 0 ~ and R-| :Rt ratios when 
the potentiometer wiper feeds a 10K load. 

In _order t p^limlnate the effects of loading^ it is jacge_conyen ient to set the po- 
t entiometers under ac tual load and monitor the wiper voltage n[th^^otentiometerout- 
put), than to calculate~a~^orrected mechanical" ratio (R*| . 

Figure 2.4-4a and b illustrate the TR-48 circuitry provided to permit setting the 
potentiometers under actual load. Relay K1 is energized -when the computer is placed 
in the pot set mode (depress PS button) and applies +10 volts reference to the Hi 
ends of all the grounded potentiometers. Note that the -wiper remains connected to 
the Pre-Patch Panel termination; thus, the viper "sees 0 the impedance of the actual 
load it is patched to in the problem. Even when the potentiometer is selected for 
monitoring via the readout system and K2 is energized, the wiper remains connected 
to its actual load. The readout system connects the -wiper to a high impedance DVM 
(10 megohms minimum) or, in the absence of the DVM, a null pot circuit. The opera- 
tor may then set the -wiper for the attenuation factor required in the problem. 

The method of setting the ungrounded potentiom eters^ is similar except the +10 volt 
reference is not automatically applied to the potentiometer high end (relay K1 is 
eliminated on ungrounded attenuators). The o perator must patch inputs to the po- 
tentiometer; this arrangement prevents possible erroneous settings depending on 
ihe~ljonfiguration in which the potentiometer is used. 

Figure 2.4-5 shows schematics and symbols for the two types of potentiometer con- 
figurations. The TR-48 potentiometer designation (i.e., number) is given within 
the circular symbol and the setting is written in close proximity to the symbol. 
In addition, the Hi and Lo ends of the ungrounded potent iometer(s) are also desig- 
nated to indicate clearly both input signal sources. 



*The high (or hi) end of a potentiometer refers to the- termination physically loca- 
ted at the top of the schematic designation of the Pre-Patch Panel. The low (or 
lo) end is the bottom termination, normally grounded except in the case of the un- 
grounded potentiometer. 
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FUNCTION SWITCHES 
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FIGURE 2.4-1 ATTENUATOR GROUP 42 283 
PATCHING AREA AND POTENTIOMETER PANEL 
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(a) INFINITE LOAD: DIAL SET AT 0,800 
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FIGURE 2.4-3. POTENTIOMETER LOADING. 
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Q. GROUNDED POTENTIOMETER CIRCUIT 
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b. UNGROUNDED POTENTIOMETER CIRCUIT 



FIGURE 2. 4- 
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POTENTIOMETER CIRCUITS , SIMPLIFIED SCHEMATIC 
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b. GROUNDED POT COMPUTER DIAGRAM 
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C. UNGROUNDED POT SCHEMATIC 
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SYMBOL 



NOTE 



k= — c 
K e 



in 



FIGURE 2.4-5. POTENTIOMETER SCHEMATIC AND COMPUTER DIAGRAM SYMBOLS 



5o OPERATIONAL AMPLIFIER 



&o General Considerations 

¥hen a high-gain d-c amplifier is used in conjunction with input and feedback 
networks to perform mathematical operations, the resulting system is generally 
referred to as an operational amplifier » The operational amplifier is the basic 
and most versatile unit in the analog computer • It can be used for inversion, 
summation, multiplication by a constant, integration, and used in conjunction 
with special networks for squaring, extracting square root, generating logarithmic 
functions, etc. 

To understand the basic concept of the operational amplifier, consider the simpli- 
fied block diagram of Figure 2.5-1 where a high-gain amplifier (gain of -A) has a 
feedback impedance Zf and an input impedance Zi n . The amplifier is designed so that 
it has three basic and essential characteristics. 

(1) The amplifier output (e Q ) is related to the summing junction voltage 
( e D ) by ' tne S^ 11 °£ amplifiers e Q = -Ae D 

(2) The input stage of the amplifier draws negligible current: lb s 0 

(3) The open loop gain of the amplifier is extremely highs A»1 (on the 
order of 3 x 10 7 at d-c). 

Using Kirchhoff's laws, the nodal current equation at the summing junction (SJ) is: 
*in 88 if + ib 

or e£ a -e D eb -eo (EQ« 2.5-1) 

_ = — + i b 

Z in Z f 

since e^ * -eo/A, and since i D « o , Equation 2.5-1 can be rewritten to obtain: 




(EQ. 2.5-2) 
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In most applications the ratio of Zf to Zi n is less than 30 and since 1/A approaches 
aero Equation 2.5-2 becomes: 



> in (EQ. 2.5-3) 



Equation 2.5-3 illustrates one of the most important considerations of the operational 
amplifier: Tfae Anput-Qtttout relationship 9f the 9PCTgUonal amplifier jg goleJLy fle- 
penflent 9R th? raj^o- of the iftreflfrftgfr ^9 th? iU^M iffP9fl«K e * 

Using Equation 2.5-3 as the basis of discussion, the following sub-paragraphs des- 
cribe the various uses of the operational amplifier. 

(1) Inversion ^ When the same value resistor is used for both the feedback and 
the input impedance, the amplifier output voltage has the same amplitude as the in- 
put voltage but is opposite in polarity. 

R f 

e o = " p e in 
R in 

In the TR-48 the value of Rf and used for the inverter is normally 100,000 ohms, 
(100K) therefore: 

100K 

e « ei n 3 -ein 

0 100K 

Thus a +10 volt input results in a -10 volt output, and the amplifier is said to 
have a gain of minus one. The accuracy of the output to input ratio depends solely 
on the accuracy of the ratio Rf/R^n. 

(2) Multiplication bv a Constants A change in the ratio of the resistors re- 
sults in multiplication by a constant. With Rf equal to 100K and R^n equal to 10K, 
for example, the amplifier output is: 

e ° = " io? eiG * ~ 10ein 

An input of plus one volt results in an output of minus ten volts. This operational 
amplifier has a gain of ten. The multiplying constant can be made smaller than one 
by using a 10K feedback resistor with a 100K input resistor. 

0 100K 2Ji 

An input of minus ten volts produces an output of plus one volt. 

(3) Summation . When multiple input resistors are used with a feedback resistor 
Rf , the basic relationship is extended to: 



U 
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FIGURE 2.5 - I . OPERATIONAL AMPLIFIER, SIMPLIFIED BLOCK 

DIAGRAM 
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DUAL DC AMPLIFIER 

b. PATCHING BLOCK 



FIGURE 2.5-2. TR-48 OPERATIONAL AMPLIFIER, SIMPLIFIED SCHEMATIC 

AND PATCHING BLOCK LAYOUT 
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FIGURE 2. 5-3. SUMMER AMPLIFIER PATCHING 



The circuit can be used to algebraically sum an indefinite number of inputs ; further- 
more, each input may be multiplied by an arbitrary constant,, 



Integration Mj&fa SSfiBSfii %9 When the feedback element Z f is a 

capacitor rather than a resistor, the summing junction current equation iss 



*1 E 2 R n dt 

Integrating this equation and assuming an initial charge on the feedback capacitor 

Looking at this another way, if is a capacitor having an operational impedance 
1/pC and Z in is a resistor, the basic operational amplifier relationship, Equation 
2 0 5-3> becomes: 



o o o 




fin - - L. '* 

pRC RC 



e_ = - -In « - L. S 0 e. n dt 



With this arrangement, the operational amplifier -will integrate (-with respect to 
time) any input voltage,, In addition to integrating, the amplifier also inverts 
the input voltage o An indefinite number of inputs may be applied to produce the 
time-integral of the sum of the input voltages o 

(5) Other Mathematical Operations o As previously indicated the operational 
amplifier has uses other than those indicated in sub-paragraph a through d« Com- 
plicated transfer functions can be simulated by using series and parallel RC net- 
works for the feedback and input impedances., The circuit performance is still 
governed by the basic relationship of Equation 2 0 5~2o For the general case where 
three-terminal networks are used, the short-circuit transfer impedance of Zf and Z in 
must be used* (The short circuit transfer impedance of a network is the ratio of 
input voltage to short-circuit output current „) The input and feedback elements 
need not be linear; therefore, almost any non-linear characteristic can be approxi- 
matedo The amplifier can also be used in conjunction with diodes and resistors to 
simulate the non-linear operations of limiting, dead-zone generation, X 2 , Log X, 
etc* 
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&o TR-4S Operational Amplifier 6o5U 

Figure 2.5-2 shows the operational amplifier patching terminations and a simplified 
schematic of the high gain d-c amplifier and summing resistor network o By placing 
a four-connector bottle plug in the patching block as shown in Figure 2.5-3a, the 
high-gain amplifier is connected to the summing resistor network as shown in Figure 
2o5~3b» The resultant operational amplifier can be used for inversion, multiplica- 
tion by a constant, and summation. The computer diagram symbol is shown in Figure 
2o5-3co On the computer diagram it is customary to show only those inputs that are 
used; the amplifier number is written inside the triangular symbol c 

Figure 2o5-4a shows the patching to provide an operational amplifier that is capable 
of integrating with respect to time. Figure 2o5-4-b illustrates a simplified sche- 
matic for an amplifier patched as an integrator,, The computer symbol for an inte- 
grator is shown in Figure 2o5-=-4co 

Figure 2.5-5 is an expanded schematic of the integrator amplifier. In addition to 
the terminations interconnected by the Tee-shaped bottle plug, certain circuits are 
brought out to the patching block for additional control of the integrators. These 
include the Operate and Reset relay coils and the operate and reset buses* Normally 
these terminations are connected as shown in Figure 2<>5-4> however, by cross patch- 
ing (operate bus to reset relay, etc.) the integrator can be used as a track and 
hold unit. 

An additional feature is the 10(3 terminations; for normal operation these two term- 
inations are jumpered by a two-connector bottle plug (Figure 2.5-4) • This bottle 
plug parallels a 9 microfarad with a 1 microfarad capacitor (Figure 2.5-5) to pro- 
duce an integrating rate of 10 volts/second for the gain-of-ten inputs and 1 volt/ 
second for the gain-of-one inputs. By removing the 1Qp bottle plug the integrating 
rate can be increased by a factor of ten. 

6 0 QUARTER-SQUARE MULTIPLIER. 

Multiplication of two variables is one of the non-linear operations necessary in a 
general purpose computer » The TR-48 Multiplier utilizes the quarter-square multipli- 
cation technique to produce a product of two variables (X and I) as illustrated by 
the following equations 

XY « - [ (X + Y) 2 - (X - Y) 2 j (EQ. 2o6-1) 



The TR-48 Quarter-Square Multiplier is basically a gated-resistor circuit applica- 
tion of the quarter-square technique. When the quarter-square multiplier is used 
as the input impedance Z± n for a high-gain d-c amplifier, the resultant circuit is 
capable of multiplication or squaring of input variables . When used as the feed- 
back impedance, the multiplier-amplifier combination is capable of division or ex- 
tracting the square root of the input variables. The quarter-square multiplier 
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b. INTEGRATOR SIMPLIFIED SCHEMATIC 

IC 



X. 



1 

1 




1 




10 




10 










= -[j*(x,+x 2 +x 3 +iox 4 +iox 5 )dt-icJ 



C. COMPUTER DIAGRAM SYMBOLS 



a. PATCHING 



4-C0NNECT0R 
BOTTLE PLUG 



NOTE 5 

IF AMPLIFIER IS NOT USED, 
BOTTLE PLUG REQUIRED TO 
PROVIDE FEEDBACK. 



FIGURE 2. 5-4. INTEGRATOR PATCHING AND DIAGRAMS 
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5. 



ONLY ONE OF TWO INTEGRATOR 
NETWORKS SHOWN. 

RELAYS K2 AND K3 COMMON 
TO BOTH NETWORKS OF 12.764. 
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ITY. (SEE PARAGRAPH 8 ). 
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FIGURE 2.5-5 



INTEGRATOR AMPLIFIER AND SIMPLIFIED SCHEMATIC 
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FIGURE 2.6-1. (X<-Y) 2 SIMPLIFIED SCHEMATIC 



employs four squaring circuits which utilize solid-state diodes and precision resis- 
tors to produce the square of a given input. Figure 2.6-1 is a simplified schematic 
typical of one of the four squaring circuits of the quarter-square multiplier. Each 
diode is reverse-biased (cut off) by the negative reference source at a potential de- 
pendent on the series resistance (R17, R23, etc.) in the reference source leg. To cause 
a given diode to conduct, the sum of the X and Y inputs (applied via weighted resis- 
tors) must attain a potential opposite in polarity and larger than the bias of the 
reference leg. By means of appropriate input resistors (R38and R37,R36and R35,etc.) 
and bias resistors (R}7,R23,etc.) the diodes may be made to conduct at different pre- 
determined values of (X + I) . (The potential at which a given diode conducts is 
known as the diode breakpoint . ) As the input (X ♦ Y) reaches the breakpoint of each 
diode, the input resistors are essentially par ailed. The gain of the operational 
amplifier changes as each diode conducts; therefore, by proper selection of compon- 
ents, the network of Figure 2.6-1 a can simulate a curve of (X + Y) 2 with straight- 
line segment s, Figure 2.6-1 b. The number of segments used determines the -accuracy of 
the approximation. 

In order to accomplish four quadrant multiplication, four squaring circuits are re- 
quired. 

Figure 2.6-2 is a simplified schematic of the quarter-square multiplier showing the 
patching terminations and the patching block area of the Pre-Patch Panel. For mul- 
tiplication or division, a four -connect or bottle plug is placed in the cross- 
hatched Q •[ area; for squaring or extracting square root the four-connector bottle 
plug is moved up to the MD area. Note in each case the symbols left exposed indicate 
the functions that may be performed. On the simplified schematic the +Y inputs are 
designated +Yj and +Yjj since they are common only when patched accordingly; these 
designations correspond to the +Y input in the upper field of the patching block 
and the +Y input in the lower cross-hatched area respectively. Resistor Rf is a 
3550 ohm precision resistor available for use as an input resistor to the amplifier 
used in the multiplier circuit. This input may be used, for example, to produce 
Z ± XY/10. 

Ac Multiplication 

Figure 2.6-3 illustrates the patching procedure for multiplication of two variables, 
X and Y. All four input (+X. -X f +Y f and -Y) are necessary even though the inputs 
flQ a<fr change 8 4g B fluriag ft profrUm rm* Note that the output voltage from the 
associated amplifier is -XY/10 (due to the inversion of the amplifier). The +X 
and -X o£ +Y and -Y inputs can be interchanged to produce +XY/10. 

h. Division 

A variable (U) may be divided by a second variable (X) by using the quarter-square 
multiplier as the feedback element of a high-gain amplifier. Figure 2.6-4 illus- 
trates the patching together with the simplified schematic for division with the 
quarter-square multiplier. 

An important point to note with this circuit is that inputs to +X must be positive 
in polarity (negative input to -X) ie # , the problem variable X, must be positive 
for division. If X is negative the +X and -X inputs should be interchanged. 
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Both +10 U/X and -10 U/X are available as outputs. 



The following restrig-UQBS m% to observed, when iig;Lng the quarter-square mvAti- 

aUfig for itolglgii» 

1o The absolute value of the divisor X must always be greater than or 
equal to the absolute value of the dividend U; i.e., U/X £ 1.0. 

2. The divisor X must be positive. 



3o The divisor X must never equal zero. 

The quarter-square multiplier may also be used to generate various other functions 
such as X 2 and«/X. The patching and associated computer diagrams for various mul- 
tiplier operations are given in Appendix 3. 

7. DIODE FUNCTION GENERATORS 

Four diode function generators (DFG's) are available for the TR-48. Three are of 
the fixed function type (X , Log X, and 1/2 Log X), and the fourth, a variable 
function type (VDFG), may be set to represent a single-valued function. All of 
the DFG f s operate on a principle similar to the quarter-square multiplier; i.e., 
diode gates are reverse biased such that the input voltage must attain certain 
levels to cause succeeding diodes to conduct. As each diode conducts the Zf/Zia 
ratio of the operational amplifier is changed. 

£• X 2 Diode Function Generator 

The TR-4.8 X 2 DFG accepts both positive and negative input voltages and generates 
either a +X 2 or a -X 2 ^ output. Figure 2.7-1 is a simplified schematic of the X 2 
generator circuit shoving the high-gain amplifiers as normally patched into the 
circuit. To generate a +X 2 function a four-connector bottle plug must be placed 
as shown in Figure 2.7-2aj this jumpers the four terminations associated vith the 
-X 2 /10 area of Figure 2.7-1$ i.e., the two upper terminations are connected to- 
gether as are the two lower terminations. For -X 2 the bottle plug is moved up 
as shown in Figure 2.7-2b. Note in each case the output function is designated 
by the exposed printing on the patching block. This feature enables the operator 
to determine the DFG output by glancing at the Pre-Patch Panel. 

An additional feature, shown on Figure 2.7-2, is the R«| input termination which 
provides an input with a gain-of-one via the amplifier connected between and 
Og. The output; shown by Equation 2.7-1 is available when an input I is applied 
to R-j. 

* j£ ± Y (EQ. 2.7-1) 

10 
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FIGURE 2.6-3. MULTIPLICATION PATCHING AND SIMPLIFIED SCHEMATIC 



A4662 




NOTES: 

1. X>U 

2. X / 0 

3. X MUST BE POSITIVE ( IF X IS 
NEGATIVE, INTERCHANGE THE 

+ X AND -X CONNECTIONS-, THE 
OUTPUTS OF AMPLIFIERS 02 AND 
03 BECOME -lOiL AND +IOii 

RESPECTIVELY ) 



FIGURE 2.6-4. DIVISION PATCHING AND SIMPLIFIED SCHEMATIC 
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FIGURE 2.7-2. X d DF6 PATCHING 



This entire quantity must be appropriately scaled to prevent amplifier overload. 

Additional functions (such as VX or X 2 /10 for limited inputs) are also obtainable 
with the X 2 DFG. The patching and computer diagrams are given in Appendix 4. 

Log X Diode Function Generator 

The TR-48 Log X DFG contains four individual log function generators all of which 
are terminated at the Pre-Patch Panel. Two of the generators accept only positive 
voltage inputs and two accept only negative voltage inputs. The output of the gen- 
erators is the scaled logarithm to the base JOji.e., 5 Log-jo [10|XI j) . 

Figure 2,7-3 illustrates the patching for both positive and negative voltage input 
log X generators. Additional patching and DFG uses (such as obtaining the antilog 
of a log input, multiplication, division, or raising an input to an unusual power; 
are given in Appendix V. 

c. . 1/2 Log X Diode Function Generator 

The 1/2 Log X DFG is similar in theory of operation, patching, and function genera- 
tion to the Log X DFG. The difference lies in the output functions the t/2 Log X 
DFG generates the log of the */% (i.e., 2.5 log 10 [10|XIJ). 

d. Variable Diode Function Generator 

Frequently in a problem, the dependence of one variable quantity (I) on another 
quantity (X) is known only in the form of experimentally obtained data. The vari- 
able Diode Function Generator (VDFG) provides a means, with a single component of 
approximating and generating functions of this type. 

The TR-48 VDFG utilizes the same technique as the Log X, 1/2 Log X, and X 2 DFG 1 s 5 
i.e., the VDFG is used to vary the Zf to Z± n relationship of the operational ampli- 
fier as a function of the input voltage to approximate a curve with straight line- 
segments. With the VDFG the operator can control only the line-segment slopes 
since the breakpoints (i.e., the beginning of each straight line segment) are fiated 
at 0, 1, 2, 3, 4, 5, 6, 7, 8, and 9 volts. 

The TR-^8 VDFG consists of two VDFG'ss a negative generator that responds to in- 



zero to +10 volts. The positive and negative VDFG 8 s can be used separately or 
combined to form a ±VDFG for both positive and negative voltage inputs. 

(1) mchjLflg and Simplified SehWttofl* Figure 2.7-4 shows the VDFG 
patching area and the simplified schematic of the overall unit; the plus and minus 
VDFG portions show only a single diode-gated resistor circuit for clarity. Figure 
2.7-5 shows the patching for using the plus and minus VDFG l s separately; this fig- 
ure also contains a simplified schematic of the plus VDFG. The minus VDFG is simi- 
lar except that the diodes are orientated in the opposite sense and are biased by 
plus reference. 
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Figure 2,7-6 shows the VDFG patching to combine the plus and minus VDFG s s for bi- 
polar inputs , Note that a four-connector bottle plug covers the 10 SEG area of the 
patching blockj one of the +IN terminations is jumpered to a -IN termination to 
combine the input leads, 

(2) Set Up Pr ocedure 0 The VDFG units are located behind the Control Panel 
as shown in Figure 2o7-7 They are individually mounted on sliding shelves thus 
facilitating access to the slope potentiometer controls* 

Each VDFG unit has eleven adjustment s; one is a parallax control and the other ten 
are slope potentiometers to adjust the output curve slope between unitary incre- 
mental inputs of X, When used in the ±VDFG combined mode the parallax pots of the 
-VDFG and +VDFG are interdependent . 

The parallax potentiometer permits the operator to set the value of I (at X = 0) 
within the range of +10 to -10 voltso The +1 slope adjustment permits the operator 
to set the initial line segment slope between 0 and 1 volt so that with 1 volt into 
the DFG, the Y output can be set to a voltage between ±2 volts above or below the 
X = 0 point (i,e,, a slope of 2 volts/volt). The remaining slope potentiometers 
(2 to 10) permit the operator the change the slope of each preceding segment by 
a voltage slope of one vgj% per yp3,t,,» 

Figure 2,7-8 is a sample output curve of a plus VDFGo This curve is used as the 
basis for the typical VDFG set-up procedure , 

The following procedure is for the set-up of a +VDFG patched as shown on Figure 
2,7-5 o The set-up adjustments of the +VDFG must be started at X = 0 and contin- 
ued in sequence to X « +10, The procedure for the -VDFG is accomplished in a 
similar manner, again starting at X = 0 and in sequence to X - -10, 

When setting up a VDFG for combined plus and minus inputs the operator again starts 
at X * 0o Once the X = 0 point is set (parallax potentiometer) the operator must 
proceed in sequence to either +10 or -10 on one VDFG, then starting again at zero, 
proceed in sequence to the limit of the remaining VDFG, 

(3) Typical YPFQ get-yp Procedure 

(&) Release the quarter-turn locking device and slide the desired 
+VDFG forward exposing the set-up controls, 

(&) Select the I = f(X) amplifier output (designated 01 on Figure 
2,7-5) for readout on the DVM (or multi-range voltmeter if the 
DVM is not available) , 

(cj Ground the +IN termination of the VDFG, Using a small-blade 
screw driver set the PARALLAX control for a DVM readout of -2 
volts. 



20 



A4666 




NOTES 

1. R, AND R 2 5000 OHMS IN ALL CASES. 

2. OUTPUTS ARE LOG TO BASE 10 



FIGURE 2.7-3 L0G. o X DFG PATCHING AND SIMPLIFIED SCHEMATICS 
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FIGURE 2.7-4 VDFG PATCHING BLOCK AND SIMPLIFIED SCHEMATIC 
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b. VFDG LOCATIONS AND ASSOCIATED PRE -PATCH PANEL AREAS 



FIGURE 2.7-7 ± VDFG UNIT LOCATION SHOWING CHASSIS IN SETUP POSITION 
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FIGURE 2.7-8. SAMPLE +VDFG OUTPUT CURVE 



(£) Apply +1 volt to the +IK termination and adjust the +1 control 
for a DVM readout of -1 0 

(&) Apply +2 volts to ♦IMo Adjust the *2 control for »1, 

(£) Continue the procedure applying X ■ +3> +4* etc <>, and adjusting 
the corresponding control for the proper D?M readout as listed 
in the table 0 

(gj For optimum accuracy repeat the set-up procedure , starting -with 
Step (&)s ana making any necessary touch up adjustments <> 

8o REPETITIVE OPERATION (REP-QP) 

The addition of the highspeed repetitive operation feature to the TR-4& provides 
a means of rapidly switching the computer integrators between the initial condition 
(or reset) and goer ate modes at controlable rates up to 5^_cjcl^_s„^ex_sec.Qnd<. The 
computer can be switched from manual to repetitive operation without requiring 
changes to the computer program.. 

Figure 2.8-1 is a simplified schematic of an integrator showing the Rep-Op capaci- 
tors and the Time Scale relay o Energizing the Time Scale relay K1, the feedback 
capacitor is changed from 10 microfarads to 0 o Q2 microfarads, thus changing the 
integrator time scaling by a factor of 500 * By terminating the time scale bus and 
relay K1 at the Pre-Patch Panels the operator can use iterative techniques in pro- 
blem solution by operating desired integrators at "real time" rates even though the 
computer is in the Rep-Op mode* 

Essentially, the Rep-Op mode, with the time scaling and rapid switching between the 
integrator reset and the operate modes, a problem is solved, reset, and then solved 
again $ this process continues repetitively at the pre -set rate» The Rep-Op circuits 
also 'provide a saw-tooth sweep voltage at the selected Rep-Op rate, suitable for 
driving an external oscilloscope 0 Thus s the problem solution may be displayed on 
an oscilloscope that is synchronized with the Rep-Op cycling rate* The sweep volt- 
age is terminated in the yellow scope area of Readout Panel 12»763 of the Pre-Patch 
Panelo (See Figure 2o3~2a<.) By connecting a bottle plug between the sweep termin- 
ation 0"O and the scope X termination of the 12°763 area, the sweep voltage is 
applied to the scope connector plug on the rear of the TR-40 (Figure 2o3~3) <> 

To place the computer in the Rep=0p mode the operator sets the desired rate on the 
REP-OP COMPUTER TIME control located on the control panelo This is a dual control; 
the outside portion is a four position switch for selecting operate times of 20 , 50 > 
1J)0 i _or__200 m illiseconds per cycle 0 The inside control (YERMIER) permits the ex- 
pansion~oi' these four settings up to 2*5 times . The control is calibrated so the 
20, 50, 100, and 200 settings are accurately obtained when the VERNIER is turned 
fully counter-clockwise , 

After setting the desired repetition rate 3 the computer is switched into the Rep- 
Op mode by depressing the R0 button on the control panel* Only those integrators 
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with a jumper bet wen the TS bus and relay K1 respond to the Rep-Op drive signals. 
9. COMPUTATIONAL ACCESSORIES 

In addition to the standard computing elements, certain devices are found useful 
in programming a problem on the computer. These items are used as required by the 
programmer to expand the computer capabilities, 

a.. Signal Comparator 40.404 

The signal comparator is an automatic switching device that consists of a comparator 
amplifier and a double-pole, double-throw relay. The amplifier compares two Input 
voltages and energizes or de-energizes the relay, depending on whether the sum of 
the input voltages is greater than zero (positive) or less than zero (negative). 
There are two separate comparator units terminated at a Dual Comparator 40*404. Br©- 
Patch Panel area. 

Figure 2.9*-1 shows the comparator patching area and a simplified schematic of one 
of the comparators. (The second comparator is identical.) The relay termination 
indicated by the negative sign Is the position of the relay when the sum of the IH-j 
and IKg inputs is negative, this is the de-energized position of .the relay. Vhen 
the sum of the IH-j and IN2 inputs is positive, the relay energizes and the wiper 
swings to the positive contact. 

(1) Comparator Switching Data, 



Relay Switching Time 

Typical * , 7 milliseconds 

Maximum allowable • 10 milliseconds 

Sensitivity 

Typical 2 millivolts 

Minimum 3 millivolts 

(2) Set Up Procedure . The problem variable that is to control the compar- 
ator switching should always be applied to the IN^ termination and the reference or 
bias voltage should be applied to the IH2 termination. The following is the set-up 
procedure for the comparator. 

(3) Apply an input to the IN* termination equal in magnitude and po- 
larity to the desired switching level (Figure 2.9-1). 

(&) Connect the wiper of a potentiometer to the IH2 termination j apply 
reference, opposite in polarity to the IH-j input, to the potent i- 
ometer hi end.'" — — 
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FIGURE 2.8-1. SIMPLIFIED INTEGRATOR SCHEMATIC SHOWING TIME SCALE CIRCUIT 
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FIGURE 2.9-1. COMPARATOR PATCHING AREA AND SIMPLIFIED SCHEMATIC 



(&) Adjust the input to IK 2 so the comparator relay is actuated as re- 
quired when IHi reaches the switching level. 

HOXE 

By applying a voltage (designated check 
voltage on Figure 2.9*-1) to one viper and 
connecting the DVM to the contact that makes 
when the relay is actuated, visual Indication 
of the relay svitching may be obtained for set 
up. 

The comparator is now set for use in the problem and should be 
patched to the appropriate signal sources. 

j£. Function Switches 

The Function Switch Group 12.766 provides five, 3-position single-pole switches 
mounted on the potentiometer panel (Figure 2.4.-5) and terminated at the Pre-Patch 
Panel area shown in Figure 2*9-2. With these switches the operator may manually 
switch computer problem circuit functions. The center position may be considered 
the srf& position. 
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FIGURE 2.9-2. FUNCTION SWITCH PATCHING AREA 



SECTION III 
BASIC PROGRAMMING 



The previous sections of this manual outlined the operation of the basic com- 
puting elements of the TR-48. An organized approach to programming a problem 
for the analog computer -will now be discussed. The sample problems -that follou 
will serve to illustrate the particular points in the programming procedure. 

1 . SYSTEM EQUATIONS 

Starting from the physical system a mathematical description can be derived using 
well known physical lavs, experimental results, and any reasonable simplifications 
that might appear desirable. The description will include a set of equations, 
coefficient values, initial conditions, and possibly some experimental data re- 
lating variables in the study. Next, it is necessary to consider the computing 
equipment available for the study, and what equation modifications, if any, are 
desirable (for instance, linearizing wherever possible to reduce the complexity 
of the study, using finite difference approximations to reduce partial differential 
equations to sets of ordinary differential equations, etc.). At the same time the 
operator should make any changes of variables or problem modifications that facili- 
tate programming. Remember that the equations most conveniently solved on the com- 
puter are ordinary differential equations. 

2. THE BOOTSTRAP METHOD 
Consider the equation 

ax + bx = c (EQ. 3.2-1) 

where a, b, and c are positive constants. 

In order to solve this constant coefficient, linear, differential equation on the 
analog computer, three devices are needed; integrators, summers and potentiometers. 

Earlier it vias shown that summation and integration with respect to time could be 
accomplished by means of operational amplifiers, and that coefficients could be 
set on attenuators. To avoid differentiation on the computer (this operation has 
inherent noise problems), the differential equation usually can be manipulated so 
that only integration is requiredo This is accomplished by re-writing the equation 
so that the highest order derivative of the dependent variable under consideration 
appears alone on the left-hand side of the equation. The rewitten equation is 

x = £ - £ x (EQ. 3.2-2) 

a a 

Note that Equation 3.2-2 indicates that if c/a and -bx/a are represented by voltages 
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and summed in an amplifier (Figure 3.2-1a), the resultant output is proportional to 
-"x (the minus sign is a result of the inversion performed by the amplifier). 

If the voltage proportional to «x is used as the input to an Integrator (Figure 
3«2-*1b) the output of the integrator will vary as x. 

The resulting output, x, may then be scaled by the coefficient b/a, multiplied by 
-1, and used as an input to the first summing amplifier (Figure 3.2-1c). 

A voltage to represent c/a is generated by multiplying reference voltage by c/lOa 

and used as the second input to the summer (Figure 3.2-1d). Thus, Figure 3.2-1d is 

a computer circuit -which, in effect, is an electronic model of the equation ail + bx = c. 

The steps folio-wed in the bootstrap procedure to obtain the model are: 

(1) Rewrite the equation to place the highest order derivative of the dependent 
variable by itself on the left side of the equation. 

(2) Assume that a voltage proportional to this variable can be formed^ con- 
sider it as the starting point in the diagram. 

(3) Integrate this voltage as many times as necessary to generate the variables 
on the right side of the Iquation. 

(4) Make use of the generated variables by. scaling them and summing them to form 
the highest derivative that was originally assumed to exist. 

(5) Provide any required forcing function and initial condition. 

Simultaneous equations may be mechanized in a similar manner as illustrated by the 
mechanization of the following set of equations (Figure 3.2-2): 

y + a-iy + a£Z + a^z = c 

(EQ. 3.2-3) 

z + b-jz + boz — bjjy - 0 

No initial condition inputs are shown to any of the integrators, so the assumption 
must be that all I.C.'s are zero* If a non-zero I.C. is needed, it can be applied 
to the proper integrator as shown in Figure 3.2-3. Bote that a sign inversion is 
present « 

One point which has not yet been mentioned is economy of equipment. In Figure 3.2-2 
the output -y from summer-1 serves no purpose other than to provide an input for the 
following integrator <> This summation might be performed by the integrating ampli- 
fier, thereby saving two amplifiers (Figure 3.2-4) . This combining of summation and 
integration may be done whenever the highest derivative does not appear anywhere 
else in the system of equations as an input and when it is not required as an out- 
put for. recording. 
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FIGURE 3.2-1 EVOLUTION OF COMPUTER DIAGRAM VIA BOOTSTRAP METHOD 
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FIGURE 3.2-3. OBTAINING INITIAL CONDITION OF y = + y Q AT t =0 
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FIGURE 3.2-4 MODIFICATION OF CIRCUIT SHOWN IN FIGURE 3.2-2 USING FEWER AMPLIFIERS 



3. SCALING 



The independent variable of the equations is represented on the computer by timej 
the dependent variables and their derivatives -with respect to time are represented 
on the computer by voltages. These signal voltages should never exceed the maxi- 
mum allowable value of ±10 volts, nor should any voltage change rapidly enough to 
exceed the frequency limitations of the computing or recording equipment . On the 
other hand, the voltages should not be so small that they are of the same order as 
possible errors in a problem solutions also, the problem solution time should not 
be so long as to result in inefficient use of the computer . Bet-ween these extremes 
there is a reasonable value for the voltage variables and the speed of solution. 
The process of scaling determines these values $ amplitude scaling to control the 
maximum value of the voltage variables, time scaling to control the speed of solu- 
tion. 

The maximum values which quantities can assume in a particular physical system can 
usually be estimated from: 

(1) A knowledge of the physical system under study. 

(2) An investigation of the system equations. Steady state relations and 
characteristic frequencies may be employed to evaluate the ranges of the dependent 
variables . 

a.. Amplitude Scaling 

For maximum accuracy, the signal voltages represexiting the dependent variables and 
their derivatives should be as large as possible without causing the equipment to 
overload. The equipment is designed to operate within the range of ±10 volts. 
This implies that the maximum expected absolute value of the problem variable must 
correspond to 10 volts on the computer. An amplitude scale factor relating volts 
on the computer to the physical units of the problem variable can then be intro- 
duced. This scale factor has the dimensions of volts per physical unit and is de- 
termined by dividing the computer reference voltage (10 volts) by the maximum ex- 
pected absolute value of the problem variable. For ease of calculation, the maximum 
expected absolute value is usually increased to the nearest convenient integer 
before calculating the scale factor 0 The table, Figure 3«3-1> shows typical scale 
factor calculations using the relationships 

10 

amplitude scale factor = i^lmated maximum of variable . 

When a range is known for a variable (for example the variable T in Figure 3.3-1 ) 
a new variable can be defined to give better resolution. By defining Z^E = T-650, 
optimum scaling for the physical quantity is obtained. 

Having determined the amplitude scale factors, the problem variables in the math- 
ematical equations are replaced by the equivalent voltages, adjustments being made 
to the coefficients throughout the equations in order to maintain equality. For 
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Problem Variable 


Max. Expected Value 
er 

Max* Expected Range 
(Physical Units) 


Scale Factor 
(Volts per Physical Unit) 


Computer 
Variable 
(Volts) 


Linear Distance, x 


0*001 inches 


— =10^ volts/inch 
.001 ' 


[10V1 


Linear Distance , s 


1 ,000 feet 


K) 1 

ITooo 53 100 volts/f 00t 


r S -j 

L 100 J 


Angular Acceleration, 9 


2-rad/sec 2 


1 n 0 
'— • = 5 volts/rad/sec 
2 


[5 © ] 


Pressure, P 


70 psi 


10 . 10 1 , . 

70 > Too = TO volts /P sl 




Temperature, T 
Let £T = T-650 


650-700 °F 
50°F 


~ - 0*2 volts/°F 
50 


[o.aff] 



Figure 3.3-1. Calculation of Amplitude Scale Factors, 
and Computer Variables. (Computer Variable*Scale Factor times Problem Variable) 



example, s and 9 "would be replaced in the appropriate equations by [s/100] and [5 9] 
respectively. The related coefficients would have to be multiplied by 100 and 1/5 to 
retain the validity of the equation. The equations are thus changed into voltage 
equations from -which a computer circuit diagram can be drawn. Summarizing the rules 
for amplitude scaling; 

(1) Determine the expected maximum absolute values and/or -working ranges 
of all variables, including any derivatives which appear in the equations. 

(2) Increase these maximum values to the nearest power of ten times, 1, 3, 
or 5 so that the scale factors can be manipulated conveniently. 

(3) Determine the scale factors by dividing 10 by the values obtained in (2). 

(4) Associate each variable with its scale factor and write the scaled volt- 
age equations from the original problem equations. 

(5) Prepare a computer diagram from the scaled voltage equations, and label 
it at appropriate points with the computer variables. 
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Time Scaling 



The independent variable o.f the computer is time and the problem solution time 
may be increased or decreased with respect to time in the physical system. Time 
scaling is carried out by defining a time scale factor p, where s 



where 



t - p t (EQo 3.3-D 

t » Time required for a solution in the physical problem. 

T = Time required foi the same solution on the computer „ 

p = Dimensionless time scale factor. 

If p is less than one, the computer solution is achieved in less time than /would 
be required in the physical system. If p is greater than one, the computer sol- 
ution takes longer than the physical problem solution. 'When (3 is equal to one, 
the computer solution is said to be a real time solution* 

The time required for a phenomenon to occur can be changed by increasing or de- 
creasing the rate at which the phenomenon takes place. This can be accomplished 
on the computer by changing the time constants in ajj, circuits by the same amount. 
Usually this requirement is met by changing the time constants of all the integra- 
tors in the problem. (If any other network associated with the problem has a 
time constant, it too must be changed to effect a change in time scale.) Changing"" 
the time constant of an integrator amounts to simply changing the "gain" of the 
integrator; increasing the input gain to decrease the time of solution, and de- 
creasing the gain to increase the time of solution. 

If the independent variable of the computer is taken as %, the integrator per- 
forms the operation 

e G (t) . - - J" e in ( T ) &% . (EQo 3-3-2) 

Since t = pt and drc = pdts 

e 0 (pt) * - jL J e in (pt) pdt = - jL- j e in (pt) dt (EQo 3*3-3) 

jr 

The integrator thus integrates with respect to time 5 t, but has a time constant 
of RC/p. The computer variable at the output of the integrator is a function 
of pt (or t). 

As an example of time scaling, consider a second order system described by 
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x + 2Cco n x + ce£ x = 0 



(EQ. 3.3-4) 



and x = 0, x = x Q , at t ■ 0 

The general analytic solution to this equation is 

x(t) - Ae <C ° nt sin (a^jvj? t + <D) (EQ. 3.3-5) 

A computer circuit to represent Equation 3.3-4 is shown in Figure 3.3-2. For the 
case where co n * 1000 and £ = 0.1, potentiometers 00, 01, 02 and the associated 
amplifiers must provide gains of 200, 1000 and 1000 respectively. These high 
gains are not compatible with maximum computer accuracy. If the circuit is pro- 
perly amplitude scaled the large gains indicate that time scaling is required 
since the rate inputs to the integrators must be reduced. 

The independent variable of Equation 3.3-4 can be changed by defining 

t = pt 

dx 1 dx 
dru " p dt 

d 2 x 1 d^ 
d T 2 p 2 dt 2 



Equation 3.3-4 becomes 



d 2 x dx ' rTrn , v 

— - + -r- + ~ x = 0 (EQ. 3.>o; 

d T 2 P dar . P 2 



whose solution is 



x( T )-Ae"^S£ sin (— N /l< 2 T + *) (EQ. 3.3-7) 

P P 

It is clear that the solution for x(<u) is simply a time scaled version of the 
original solution for x(t) with the natural frequency (aO and decay time (Ccoq) 
reduced by a factor p. The amplitudes and initial conditions are not affected. 

An appropriate time scale factor for this example is 100 which results in a com- 
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-10 




FIGURE 3.3-2 COMPUTER DIAGRAM FOR EQUATION 3.3-4 




+ X 



10/3 



FIGURE 3.3-3. TIME SCALED DIAGRAM FOR EQUATION .3.3-4. 
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COMPUTER DIAGRAM 
-I0<X<0 OR 0<X< + I0 
-I0< U<+I0 BUT ju| < |X| 

DIVISION CIRCUIT RESTRICTIONS 

1. #< 1.0 

2. X MUST BE POSITIVE, (IF X IS NEGATIVE, 
CONNECTIONS TO +X AND -X SHOULD 
INTERCHANGED) 

3. X # 0 



a m-i 



puter frequency of (co n /p) of 10 rad/sec and a damping time constant (Coon/p) of one 
second. This change in time scale is obtained as shown in Figure 3.3-3 by dividing 
all inputs to the integrators by p. Note, in addition, that potentiometer settings 
less than or equal to one are obtained by using gains of ten on the integrators. 

Let us check through the circuit of Figure 3.3-3. 

The input to amplifier 01 is 

6i1 " iqp 

The output of amplifier 01 is 

e o1 =-10 J dLtd] = 5 xdt = *(,) 

1 Up 

The input to amplifier 00 is 

The output of amplifier 00 is therefore 

e o0 = J [pi *(t) + r x(t)] d T = H- — ) * - - — 

Thus with p - 100, the rates to the integrators are decreased by 100, the computer 
solution is slowed down by 100, and the time scale is effectively 100 times its 
previous value. 

Since the amplitudes present in the computer solution are not changed by time 
scaling, it is reasonable to use the same notation for the variables in the time 
scaled circuit as is used in circuits without it. It is important, however, to 
remember that all computing components, including output equipment such as strip 
chart recorders and oscilloscopes, are operating in machine time (or %) and the 
time axis of all recordings must eventually be labeled in terms of t in order to 
interpret the output data. If for example, a particular point in the time scaled 
computer solution is labeled as ten seconds (in t) } then the same point is labeled 
as 10/p seconds in t. 

The procedure for time scaling may be summarized as follows? 

(1) Consider in the early stages of the problem whether or not a time . 
scale change is indicated. If possible, determine the natural 
frequencies and exponential time constants of the physical system. 

(2) Perform amplitude scaling and write the scaled voltage equations as 
described earlier without regard to time scaling. 
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(3) If the input gains to the integrators are very high, or very low, 
include in aJJ, the integrator inputs an attenuator factor of 

(4-) Choose an appropriate value of p from steps 1 or 3, and calculate 

the necessary potentiometer settings. If possible, select (3 so that 
standard amplifier gains and reasonable potentiometer settings (>.Q1) 
are used. 

(5) Determine the recorder speed settings or scale factors for output 
devices so that time axes read appropriately. 

Uo COMPUTER CIRCUIT DIAGRAM 

A computer circuit diagram can be produced directly from the scaled voltage equa- 
tions. It -will follow the general data flow of a mathematical block diagram but 
•will include detailed information about the computer components such as, compon- 
ent assignment number, amplifier and integrator gains, potentiometer settings 
(in terms of the system parameters) and function switch terminations, etc. The 
outputs of major components should be labeled with the computer variable appearing 
at that point in the diagram. Standard symbols should be used in preparing the 
diagram such as the recommended set given in Appendix 1. The diagram should be 
drawn so that it can be easily reproduced. 

Each computing component on the diagram is assigned a number corresponding to the 
location of the component in the computer. The assigned numbers appear on the 
amplifier assignment sheets and the potentiometer and function switch assignment 
sheets. The amplifier assignment sheet should include problem identification, 
amplifier number, function of the amplifier (summer, inverter, high gain, or 
integrator), the variable represented by the output of the amplifier, calculated 
values for the static check voltages and provision for recording the measured check 
voltages and any comments or notes. The potentiometer assignment sheet should in- 
clude problem identification, potentiometer, number, parameter description, poten- 
tiometer settings for both static Qheojc and the problem runs. Function switch 
assignment sheets are prepared in essentially the same form as potentiometer as- 
signment sheets. 

5. PROBLEM CHECK PROCEDURES 

When the programming and scaling operations are completed and the computer cir- 
cuit diagram is drawn, the methods to be used in checking the computer program 
and subsequent patching should be considered. 

Amplifier and potentiometer assignment sheets are' prepared and a detailed sequence 
of steps to be followed after insertion of the patch panel should be drawn up. 

The first phase of a problem check procedure is the static check. The static 
check consists of two parts. Part one, the program check is a pencil and paper 
check to determine if the scaled computer program (circuit diagram) truly re- 
presents the original equations. This part should be completed before patching 



31 



the problem. Part two, the .patching check, is the actual measurement of com- 
puter input and output voltage levels to ascertain that the interconnections 
(made by patching) of computer components are correct, and to check the static 
operation of inverters, summers, high gain amplifiers, potentiometers, multipliers, 
X DFG's, etc. The static check does not test the operation of the integrators. 

To perform the first part of the static check, assume convenient values for the 
variables, their n-1 derivatives, and all parameters. The variables and their 
derivatives are represented by scaled initial condition voltages at each integrator 
and need not have physical significance. The values selected for these initial 
conditions should be chosen to produce reasonable outputs (one to ten volts) from 
all computing components. The selected scaled voltages are then written down on 
the computer diagram next to the outputs of the appropriate integrators. These 
voltages are then used as a starting point for computing the voltage levels of all 
input and output terminations of the computing components on the circuit diagram. 
These computations are made by performing the various operations of multiplication 
by a constant, sign inversion, summations, etc., at the computing components on 
the diagram until the input to an integrator is reached. At the input to an inte- 
grator the calculation stops since the integrator is not operational in the I.C. 
mode. This procedure is repeated for all variables until the complete set of volt- 
age levels has been calculated. 

The assumed check values are then substituted into the original unsealed mathemati- 
cal equations and the value of the highest derivative is computed. This computed 
value for the highest (nth) derivative, properly scaled, (includes both the voltage 
scale factor^ and the time scale factor R) should check the voltage level calcu- 
lated for that point on the diagram. In order to check the inputs to the other 
integrators, a calculation is made on the basis that the input to the integrator 
should be minus the derivative of the output (divided by R, the time scale factor). 
Example: If the output of an integrator was [-1000 x] the input should be +1000x/p. 

If the assumed value of x = 0 o 05 ft/sec^ and R = 10 this results in an input volt- 
age of +1000(.05)/10 or +5 volts . The +5 volts should check the calculated input 
voltage level on the diagram. In some cases the assumed values for the variables, 
n-1 derivatives and parameters produce low voltage levels or zero outputs. In 
these cases it is necessary to change some of the assumed values and repeat the 
calculations. In the case of coupled equations a static check should be made of 
the entire system with one set of assumed values. When this portion 'of the static 
check is complete, the computer program should be patched „ 

Patching is best accomplished by two people, one person calling off from the com- 
puter diagram, the connections to be made, and the second person making the con- 
nections on the patch panel. A colored line is then drawn through the corres- , 
ponding line on a copy of the computer diagram indicating that a connection has 
been made. After patching is complete the two persons should exchange places, 
and the patched program read off from the patch panel by one person and checked 
against the diagram by the other. A systematic procedure that is recommended is 
to check out all of the patched inputs to computing components against a copy of 
the computer diagram. 
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The second part of the static check is to set the potentiometers and initial con- 
ditions on the computer to the static check values. Since all computing compon- 
ents except integrators are operational in the I.C. mode, their outputs can be 
read out and checked against the calculated values. In order to check inputs to 
integrators it is necessary to use an uncommitted amplifier. The inputs to the 
integrators are temporarily patched into this amplifier to allow read-out of 
the effective input voltage to the integrator. Since the input to an integrator 
can exceed ten volts it is often necessary to use a check amplifier -with a gain 
of -1/1Q. The use of a check amplifier is described in the example problems. 

Once the static check is completed one or more dynamic checks are made depending 
on the information available. In linear systems a simplified form of the equa- 
tion can be checked by reducing the damping to zero and checking undamped natural 
frequencies. Many times the response of the system to a step input can be cal- 
culated and the computer program modified to accept a step input as a check. 
In non-linear systems, a linear approximation can often be implemented by adjust- 
ment of the parameters and a check obtained using results from a precalculated 
linear analysis. Numerical checks from digital computer results can also be used. 
If experimental data for the model is available, a check can be made by direct 
comparison. 

6. PROGRAMMING A LINEAR PROBLEM 

To demonstrate programming, amplitude scaling, and time scaling, a computer pro- 
gram will be developed to investigate the effect of shock-absorber damping coef- 
ficient upon the transient motion of an automobile chassis. 

a» Problem Description 

A simplified representation of one wheel of an automobile suspension system to- 
gether with one-fourth the vehicle chassis, is shown in Figure 3.6-1 . 

The differential equations of motion are derived by equating the forces acting 
upon the masses to zero to establish dynamic equilibrium. 

Thus for M-js 




Z 



(EQ. 3.6-1) 



And for M 2 : 




+ K 1 (x 2 - x^ + K 2 (x 2 - x 3 ) u 0 (EQ. 3c6-2) 



Inhere 



M. 



i 



g 



= 25 



slugs K-| = 1000 lb/ft 
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FIGURE 3.6-1 SIMPLIFIED REPRESENTATION OF A SINGLE AUTOMOBILE -WHEEL SUSPENSION SYSTEM 
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FIGURE 3.6-2 MATHEMATICAL BLOCK DIAGRAM FOR EQUATIONS 3.7-3 AND 3.7 



M 2 = — = 2 slugs K 2 = 5000 lb/ft 



= x(t) = step inputj 5 inches in amplitude (car riding up on curb) 

D = variable damping coefficient | 20, 50, and 100 Ib-sec/ft 

For scaling -we shall assume that the displacements x«j and x 2 "will not exceed twice 
the input function and that the velocities x-| and x 2 will not exceed 5 feet per 
second and 50 feet per second respectively » 

Preliminary Considerations 

The mathematical block diagram is the first step in programming » It is a trans- 
lation of the system equations and shows the basic approach to the -whole problem* 
From this diagram one can determine approximate equipment requirements! it also 
serves as a guide to the final form of the scaled equations , the method of gen- 
erating the various functions required, and the layout of the final circuit dia- 
gram. 

The block diagram shown in Figure 3.6-2 was obtained from Equations 3° 6-3 and 
3.6-4 using the bootstrap method* The initial equipment requirements are four 
integrators, four inverters and 10 potentiometers (exclusive of IC potentiometers) 

X 1 = " Ej" *1 + TEf *2 " x 1 + 17 x 2 + h (EQ- 3 .6-3) 

An examination of the natural frequencies and damping time constants of the un- 
coupled Equations, 3.6-5 and 3. 6-6, will give a good indication as to whether 
time scaling will be required . 

Dx 1 K 1 x 1 

* 2 + W + MSTT' 2 " 0 (EQ " 3,6_6) 

Inhere the undamped natural frequencies ares 



co^ = J^- - »JZ0 = 6,325 rad/sec = 1 cps 
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°° n 2 = V ~^ = *y35UcT » 54.77 rad/sec - 9 cps . 



Using the maximum value of D, the damping time constants ares 

2M 1 2M 2 

= 0.50 sees — — = 0o04 sees 

D D 

Therefore, for the uncoupled systems, the frequencies involved will be of the 
order 1-9 cps, and the time required for the solution to reach a steady state 
■will be in the order of 5 time constants or 0.2 to 2.5 seconds* Frequencies in 
the order of 0.1 to 1 cycle/second are more desirable, therefore, a time scale 
change with 10<p<100 is indicated. 

c,. Scaling 

The amplitude scale factors are obtained by dividing 10 volts by the maximum 
absolute value of the problem variables. The computer variables are obtained 
by multiplying the problem variables by their scale factors. In this example 
the maximum values are given, and the computer variables are calculaied in the 
table, Figure 3.6-3 « Since the accelerations, x-j and x 2 are not explicitly re- 
quired in this study, these terms can be formed at the inputs to integrators and 
can be scaled at the level of x-j and x 2 respectively. 



ftrobleui 


Maximum 




Computer 


Variable 


Value 


Scale Factor 


Variable 


x 1 


10/12 feet 


120/10 10 volts/ft 


[1Q*1] 


x 2 


10/12 feet 


120/1 0> 10 volts/ft 


[10x 2 J 


*1 


5 ft/sec 


10/5> 2 volts/ft/sec 




*2 


50 ft/sec 


10/50> 0e2 volts/ft/sec 


[0 a 2x 2 j 


• • 

X-j 


Use same S.F„ as for x-j 


2 

2 volts/ft/sec 




». 

*2 


Use same S.F. as for x 2 


0.2 volts/ft/sec 2 


[0.2x2] 


x 3 


constant 







Figure 3«6-3 0 Calculating Scale Factors 
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Solving the original equations for the highest derivatives and introducing the 
scaled variables yields the scaled voltage equations-. 

TJ i D 1 ^1 1 ^1 "I 

[0.2x 2 ] = -0.2 — ^ [0.2x 2 J +0.2 - [2^] -0.2 — ~ [1Qx 2 ] +0.2 J£ ^ [10^] 

-0.2 JL [10x?J +0.2 i. [-jo] - .Z. r .--r- r ? (EQ. 3.6-8) 

' M 2 10 * M 2 10 



Mote that square brackets are placed around the computer variables. The computer 
variables are voltages which have been scaled to be as large as possible but less 
than 10 volts . These equations contain constant coefficients •which will appear in 
the computer circuit as, attenuator settings (that cannot be greater than one) and 
amplifier gains. The maximum value of these gains can be determined after first 
simplifying Equations 3.6-7 and 3.6-8. 




(EQ. 3.6-10) 



These gains must be obtained by a combination of potentiometer settings and ampli- 
fier gains o 

Tne circuit to generate [2^] is shown in Figure 3.6-4-a where standard gains of 
one and ten are used where possible and the attenuator settings are less than one. 
The nigh input gains required suggests that the input rate is high; it is an indi- 
cation that time scaling is required. The rate can be decreased by a factor of 10 
to slow tne problem down. The circuit with p = 10 is shown in Figure 3.6-4-b. 
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Including p , the final scaled equations are : 

L— J ■ -KjJJp)[2x t ] + 10(— )[0.2x 2 ] -1 (^pltlOx^ +1 i (5M^p)[l0x 2 ] (EQ. 3.6-11) 

0 2x K K 

[ r J]= " 10( W [0 - 2i2] +1 ( W C2il3 " 1( # )t,0X2] +1( <% )Cl0x i ] 

-10(-= )[10x J +10( J )[1Q] (EQ. 3.6-12) 

500M2p 2 J V 500M2p L J 



In Equations 3.6-11 and 3.6-12, the terms in square brackets represent computer 
variables; the terms in curved brackets represent attenuator settings.; the un- 
bracketed coefficients represent amplifier gains. From these equations the com- 
puter circuit diagram can be drawn directly. 

Computer Circuit Diagram 

Taking Equation 3.6-11 as an example: 

(1) Integrate the highest derivative to obtain the lover derivatives and 
the variable itself; Figure 3.6-5. 

(2) Similarly generate -0.2x 2 and 10x 2 . 

(3) Use the scaled equations to provide the inputs to the first integrator, 
e.g., the first term in Equation 3.6-11 says to take -2x^ through an attenuator set 
at D/M-jP into a £ain of 1 on the integrator; the second term says to take +0.2x 2 
through an attenuator set at D/M-jP into a gain 10 on the same integrator. The other 
terms in the x-j equation are supplied in a similar manner. The complete computer 
diagram is shown in Figure 3.6-6. Note that integrator 06 is generating +[0.2x 2 ], 
rather than -[0.2x 2 ]. This circuit arrangement uses two less amplifiers than the 
circuit that -would follow directly from the mathematical block diagram. Also note 
that numbers have been assigned to all components and all amplifier outputs are 
labeled with the appropriate computer variable. 

e. Check Procedures and Set Up Sheets 

Once the scaled computer diagram is completed, a static check is calculated and 
the equipment assignment sheets are prepared. A definite sequence of steps must 
be followed in preparing a static check so that the operator can check his work 
up to this point. 

(1) Calculate all potentiometer settings, using convenient values for 
the parameters, say those which will be used in the first run. For this example, 
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+ IOX 2 

a. DIRECT MECHANIZATION OF EQUATION 




d(2Xi) 
d T 



b. TIME SCALED CIRCUIT 



FIGURE 3.6-4 DIRECT AND TIME SCALED MECHANIZATION OF EQUATION 3.6-9 



-IOX, 




NOTE 

A GAIN OF 5 IS REQUIRED THROUGH 
THE INTEGRATOR TO OBTAIN IOX, 
FROM 2X, AND THE FACTOR OF £ 
MUST BE INCLUDED TO TIME SCALE 



FIGURE 3.6-5 GENERATION OF IOX, FROM 2X, 
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FIGURE 3.6-6 SCALED COMPUTER DIAGRAM FOR THE AUTOMOBILE SUSPENSION SYSTEM 



choose the parameters as given and D = 100 Ibs-sec/ft. List the parameter and 
static check setting of each potentiometer on a potentiometer set up sheet as 
sho-wn in Figure 3.6-7 . 

(2) Assume convenient values for the inputs and initial conditions. 
These values need not have physical significance. For this example choose: 

x^ = +1 ft x-j = +5 ft/sec x^ = 3 inches 

x 2 = -1 ft x 2 = -10 ft/sec 

On a copy of the computer diagram^ "establish" these values as initial conditions 
or inputs, see Figure 3.6-8. If additional potentiometers are necessary, list 
them on the potentiometer set up sheet. 

Note that integrators that normally have a zero initial condition for problem 
runs are no-w supplied -with initial condition inputs. ffhese inputs must be re- 
moved -when the static check is comp leted. The reference voltage to potentiometer 
15 and to the initial condition inputs of amplifiers 02, 03, and 07 -will be supplied 
by the reference terminations next to potentiometers P15-19° The minus reference 
from this source will be designated -6Tj the positive reference, +ST . "When the 
static test is completed, all leads to this reference block -will be removed. 

(3) On the basis of the initial conditions and inputs established, deter- 
mine the voltage at the output of each computing component as showi in Figure 3.6-8. 
The output of A03 is +10 volts s the output of A00 is -10 volts, the output of P02 

is -8 volts, and so on. This checks that reasonable outputs are being obtained 
from all components for the values chosen. If an overload or zero voltage results 
from anv component . choose different initial conditions. Also list the output 
voltages of all the amplifiers on the amplifier set up sheet as shown in Figure 
3.6-9. 

From the diagram, determine the voltages representing the highest derivatives, e.g., 
in Figure 3.6-8, the voltage representing x 2 is 

-0.2xo 

- = -5(1) -10(1) -1.25(10) -5(10) -10(1) -1(10) = -97.5 volts 



In a similar fashion; 
+2x 1 

= -28 volts (integrator 02) 

H 

-10x 1 

■ = -5 volts (integrator 03) 
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- = -10 volts (integrator 07) 

P 

These voltages will appear at the inputs to the integrators and can be measured 
by making use of a check amplifier as explained later. List these voltages next 
to the appropriate integrator under the calculated check point column on the 
amplifier set up sheet. 

(4) Substitute the assumed initial condition values of the variables into 
the original system equations, calculate the value of the highest derivatives and 
then the scaled computer vol cages to represent these derivatives. These voltages 
should check with the voltages determined in step three. Substituting into Equa- 
tions 3.6-3 and 3*6-4: 

' x 'i - If (5) + W ( " 10) - IP (1) + IP ■ - uo " /sec2 

x 2 = 4875 ft/sec 2 

Therefore the input to integrator 02 should be 

. &=mi m volts 

p 10 
the input to integrator 06 must be: 

±£g . ~Q>2(4B75) . _ 97 . 5 volts 
p 10 

This check verifies that the final program is a model of the original equations. 

(5) The final step in the static check is to actually measure the volt- 
age outputs of the computing components. Any discrepancy between calculated and 
measured outputs must be due to incorrect potentiometer settings, faulty patching, 
or to a component malfunction. 

After the problem is patched from the computer diagram, the static test inputs are 
introduced. The computer is placed in the Pot Set mode and the potentiometers are 
set to the coefficients listed for the static test (ST) case on the potentiometer 
set up sheet. The computer is placed in the RS mode and the voltage output of 
each amplifier is measured and recorded on the amplifier set up sheet. These 
measured values should check with the values calculated in step three. 
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TR-48 POTENTIOMETER ASSIGNMENT SHEET 



DATE PROBLEM 



POT 
NO. 


PARAMETER 
DESCR IPTION 


SETTING 
STATIC 
CHECK 


SETTING 
RUN NO.I 


SETTING 
RUN NO. 


SETTING 
RUN NO. 


NOTES 


POT 

no; 


00 


5//3 


.5000 


.5000 




• 




00 


01 


D/M./3 


.4000 


.4000 








01 


02 


K./5M, jS 

_ . , . ; 


.8000 


.8000 








02 


03 


j 

K,/50M 2 /3 


1.000 


i.ooo ; 






03 


04 




! i 

1 






04 


05 


D/M./S 


.4000 


, , 

.4000 








05 


06 


5//S 


.5000 


.5000 








06 


07 


D/IOM 2 £ 


.5000 ' 


.5000 








07 


08 


D/IOM 2 /3 


.5000 


.5000 








08 


09 














09 


10 


K,/50M 2 /3 


1.000 


1.000 








10 


'I 1 


K 2 /500M 2 /3 


.5000 


.5000 








1 1 


12 


K,/5M iy S 


.8000 


.8000 








12 


13 


K 2 X 3 /500M 2 /9 


.1250 


.2083 






X 3 = 3 " FOR ST 


13 


14 








j _ j 






14 


15 


0.2 X 20 /I0 


.2000 


REMOVE 






STATIC TEST I . C. 


15 


16 














16 


17 














17 


18 














18 


19 














19 


20 












- 


20 


21 














21 


22 


— 












22 


23 














23 


24 














24 


25 














25 


26 














26 


27 














27 


28 














28 


29 














29 



FIGURE 3. 6—7 TR-48 POTENTIOMETER ASSIGNMENT SHEET 
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NOTES: 

1. ALL PATCH CORDS TO ST REFERENCE 
ARE REMOVED WHEN CHECK IS 
COMPLETED . 

2. VOLTAGES IN RECTANGLES ARE 
STATIC CHECK VOLTAGES. 



FIGURE 3.6-8 STATIC CHECK PROGRAM 
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TR-48 AMPLIFIER ASSIGNMENT SHEET 
DATE PROBLEM 



AMP 
NO. 


FB 




STATIC CHECK 


NOTES 


FUNCTION, AND/OR | 
VARIABLE 


CALCULATED 


MEASURED 




CHECK PT. 


OUTPUT 


CHECK PT. 


OUTPUT 


00 


S 


-IOX, 




- 10 




- 10 




01 


IOK 


CHECK AMPL. (-1/10 ) 












02 


INT 


-2X, 


- 28 


- 10 


+ 2.8* 


- 10 


* CHECK AMPLIFIER 
GAIN - 1/10 


03 


INT 


+ 1 0 X | 


- 5 


+ 10 


+ 5 


+ 10 




0 4 
















05 


S 


4-10 X 2 




- 10 


T i 
\ -10 1 


06 


INT 


+ 0.2 X 2 


-37.5 


- 2 


l — r — j 

+ 075 |-2 

4 __j . 


O^ 


INT 


— 10 X 2 


- 10 


+ 10 


+ 10 | +10 | 


oa 


i 

! 






1 i 


09 










s 

111 

i i i 


10 










i 


; 




II 










I 


■ 

i 


12 










| 




13 










1 ! 

1 \ _ 


i 

; 1 


14 










i 
I 

1 




15 












: ' 1 


16 












17 














18 










' y . ' 

I ! 'A 


19 










\ '■> i 

- __ . _ f .„ - i. .... 


20 




— ' : ■ 






f 

i 


! ! 

! 1 
] \ 


21 










f- 

1 

i 


22 










1 — — — f - - 

j 


23 










i : 

i l 



FIGURE 3.6-9. TR-48 AMPLIFIER SETUP SHEET 



The check points of the integrator can be quickly checked by making use of a 
"check amplifier". For example, remove the SJ to SJ connection for integrator 
02 and connect the summing junction of the four input resistors to the base of 
amplifier 01 (-which, in this case, is not being used in the problem). Since 
amplifier 01 has a 10K feedback resistor, its gain is -1/10 and its output should 
measure +2.8 volts. The circuit is shown in Figure 3.6-10. 

If all the voltages on the patched computer program check against the calculated 
values, then the patching, pot settings, and amplifier gains are verified. 

If the check point voltage of an integrator is less than ten volts, the gain of 
the check amplifier is increased to -1 „ The measured check point values and the 
gain of the check amplifier are listed on the amplifier set-up sheet. 

7. A NOH-LIHEAR PROBLEM 

The next example is supplied through the courtesy of the Engineering Department 
of E.I. DuPont de Nemours and Company, Wilmington, Delaware • Again the primary 
concern is to get the problem onto the computer. 

&o Problem Description 

A surge tank must accommodate fluctuations in flow of effluent from a plant* The 
outlet line from the tank is terminated by a weir, and the flow goes to a river. 
The system can be represented as shown in Figure 3°7-1o 

The tank is in service and operates satisfactorily with two pumps feeding into 
the tank. Expansion of the process calls for the addition of another pump. A 
steady state calculation shows that the tank will not overflow with all three 
pumps operating. The question is raised whether the tank might overflow before 
the mass of water in the outlet line could accelerate to the final discharge 
rate*, If this occurs, additional outflow or tank capacity must be provided. 
Operating considerations indicate that all pumps have the same capacity and that 
the third pump will not start until well after the first two have started. 

The system is described by three equations? 



Material Balance in Tank; 



(ft 3 /sec) 




(EQ 0 3.7-1) 



Force Balance? (lb) 




(EQ. 3o7-2) 



Empirical Formulas (ft) 

For Head Loss Due to Pipe Friction 



H p = 0.0185 Q L + 0.000036 Q L (EQ. 3*7-3) 



4-0 



N = number of pumps operating 

= flow from one pump, ft^/seq 
Q L ■ outflow from line, ft^/sec 
Hrp = height of liquid in tank, above weir, ft 
Hp « friction head loss in pipe, ft 
t = time, sec 

Aj, ■ cross sectional area of tank, ft 2 

Ajj == cross sectional area of line, ft 2 

L = length of line, ft 

P » density of effluent, lb/ft 3 

g = gravitational constant, ft/sec 2 
The system, parameters are: 

Ap = 470 ft 2 Q A = 90 ft 3 /sec 

A L = 12.6 ft 2 N = 1, 2, 3 

L = 665 ft Bj < 9 ft, to avoid overflow 

g = 32.2 Ft/sec 2 

Design a computer program to study this system, switch in the pumps with function 
switches. Does the present tank have sufficient capacity? If not, how high 
must the tank be? Are there any other possible solutions aside from increasing 
the height of the tank? 

Preliminary Considerations 

The three system equations can be reduced to two by substituting Equation 3.7-3 
into Equation 3.7-2. Rewriting the resulting equations so that the highest 
derivative appears alone on the left-hand side: 

\ = (uq a - q l ) (EQ. 3.7-4) 

Q L = ftffe (H T -0.0185 -0.000036 q£) ( E( i- 3.7-5) 

Li 
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_4V 



- 8V- 



-0.8 V 



-8V 



! I00K 
— VW- 



10 



IOOK 
-AAAr- 



I0K 

-AAAr 



100 K 
■AAAr 



MAKE THIS CONNECTION 



INPUT NETWORK 
OF A02 



2.8V 




CHECK AMP. 



FIGURE 3.6-10 CHECK AMPLIFIER PATCHING TO MONITOR INTEGRATOR OR CHECK POINT, TYPICAL 
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FIGURE 3.7-1. 



SURGE TANK, SIMPLIFIED SYSTEM DIAGRAM. 
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-10 









Qa 




2Qa 




3Qa 


I0A T 


I0A T 


IOA T 




A 

1 a o 1 





— NQ* 



O.OI85gA L 



y _ 0.000036gA L 



I 

A T 



+ Qi 



-i 



-/dt 



+ H T 



9A, 



-/dt 



FIGURE 3.7-2 MATHEMATICAL BLOCK DIAGRAM 



The mathematical block diagram is shown in Figure 3.7-2. One possible arrange- 
ment for switching the pumps is shown. Note that Q^ is a constant and need not 
be scaled o 

Equation 3.7-5 can be used to determine the steady state level of the head (HQ^ = 
QjJ with the following results: 



N = 1 


Hj = 


1.96 ft 


N = 2 




4*5 ft 


N = 3 


% = 


7.62 ft 



The computer results should confirm these steady state levels. 
£. Scaling 

The maximum value of the inflow is given as 270 ft^/sec (3 pumps on). The maxi- 
mum value of Qjq is not known; it may be greater than the inflow due to surges in 
the tank. The maximum value of Ql will be estimated at about twice that of the 
maximum inflow, say 500 ft^/sec. This value may be changed later if the ampli- 
fier voltage representing Ql is extremely high or low. 

The maximum allowable value of H^, before tank overflow, is 9 feet. If % is 
scaled at this level and the tank overflows due to surges, the amplifier pro- 
ducing Hx will overload • It would then be impossible to determine the amount 
of additional tank capacity required. Id order to avoid this difficulty, assume 
that extra capacity already exists, i.e., the tank is about twice as high, say 
20 feet. Estimates of % and §l are not required since these variables are of 
no interest and will appear only as inputs to integrators; they may be therefore, 
scaled at the same level as % and Ql. The amplitude scale factors and the com- 
puter variables are shown in Figure 3*7-3* 

The scaled equations are obtained by substituting the computer variables into 
Equations 3*7-4 and 3*7-5 » 



Variable 


Est* Maxo Value 


Scale Factor 


Completer Variable 


Ql 


500 ft^/sec 


10 
500 


[0 o 02 Q L ] 




20 ft 


10 
20 


[0.5 %] 


%, 


Scale to Level 




[0-02 Qj 


o 


Scale to Level 


of % 


[0o5 %] 



Figure 3*7-3* Scale Factors 
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[0.5HT] = m A - §^ [0.02Q L j = 2^|10| -mg [0.02QJ (EQ. 3.7-6) 

("Where the ten in braces refers to reference voltage which is used to generate 
the constant term NQa/At.) 



rn nV, 1 °- 02 g A L rn ™ 1 °-W°-W)&L r n ma , 0.02(0.000036)^ ; 
[0.02Q L J = — — [0.5%] - — [0.02QJ - [ — ) (Q _ 02)L [0.02Q L ] 

[0.02QJ = ^ [0.5* j - [coafcj - 2=22S5Sk [£^M£^ 

(EQ. 3.7-7) 

(The square term for Ql is arranged so that the scale factor of 1/10 inherent in 
multiplication is considered.) 

The low value of the coefficients in Equations 3.7-6 and 3.7-7 indicates a long 
solution time, therefore, time scaling is desirable. Each of the terms on the 
right-hand side of these equations appears as an input to an integrator. Time 
scaling is accomplished by dividing all integrator inputs by an appropriate factor 
p, where p is selected to give reasonable pot settings. In this case let p * 0.1. 

Computer Diagram 

The computer diagram is shown in Figure 3.7-4-. The function switches make it 
possible to switch in the pumps, one at a time, in case a particular timing se- 
quence will give a more satisfactory response. Four gains of one are needed for 
integrator 02. This can be accomplished by borrowing a resistor from an unused 
amplifier and connecting it to the summing junction of the integrator. Since Q2 
can never be negative, only one squaring card of the X 2 DFG is used. Integrator 
15 is used to provide a time-base for plotting % as a function of time on an 
X-I plotter. 

e. Checking the Computer Program 

Following the steps used in the previous example; 

(1) Calculate all potentiometer settings and list them on the potentiom- 
eter set-up sheet. 

(2) Assume convenient values for the variables. 

lip = 20 ft 

Q L = 500 ft 3 /sec 

= 90 ft^/sec (one pump on) 
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1000/3 



FIGURE 3.7-4 COMPUTER DIAGRAM; SURGE TANK 



Establish these values on a copy of the computer diagram. On the basis of the 
inputs and initial conditions established, determine the voltage at the output 
of each computing component, and the voltages representing the higher derivative, 
The voltages representing Hj and at the input to integrator 02 and 03 respec- 
tively are: 

—Si = -10( .096)(1) +10( .532)0) = +4.36 volts 
P 

+0.02Q r 

2k = +2.440) -1.130) -1.10(1) = +0.21 volts 

P 

(3) Substitute the assumed values of the variables into Equations 3.7-4 
and 3.7-5, calculate the values of the highest derivatives. 

= ^ (90 - 500) = -0.8723 ft/sec 
q l = ^32.2)(12.6) [20 - 0.0185(500) - 0.000036(500) 2 ] = 1.067 ft 3 /sec 2 

therefore 

= ^(-0.8723) m ^ volts 

p 0.1 

+Q,Q2 ^ - lO^Lm = +0 .213 volts 
p ' 0.1 

•which checks the computer diagram. 

(4) Check that the outputs of the amplifiers are labeled correctly by 
substituting the assumed initial values into the expression for the amplifier 
output. Prepare amplifier set-up sheets. 

(5) Place the computer in the RS mode and verify that the calculated and 
measured outputs agree. This checks the patching, pot settings, etc. 

(6) Remove all static test reference leads, check that the integrators 
are reconnected properly and prepare for run one. 

The steady state values of % for N = 1, 2, 3 were determined earlier; the com- 
puter should confirm these levels. Several other check cases can also be deter- 
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mined. For example, assume the system is in equilibrium with one pump operating $ 
then pump two is switched on. Maat -would be the instantaneous rate of change of 
the head (ft/sec) in the tank? 

Consider the following possible solutions ; increase the height of the tank, increase 
the tank diameter, increase the pipeline diameter, drive the third pump with a two- 
speed motor, etc. (The addition of a third pump identical to the first two may be 
due to plant standardization and the pump capacity may be arbitrary.) Find the 
largest capacity that can be added and still prevent overflow. 
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TR-48 AMPLIFIER ASSIGNMENT SHEET 
DATE PROBLEM 



AMP 

NO. 


FB 




STATIC CHECK 


NOTES 


FUNCTION, AND/OR . 
VARIABLE 


CALCULATED 


MEASURED 




CHECK PT. 


OUTPUT 


CHECK PT. 


OUTPUT 


00 
















01 










i 






02 


INT 


+ 0.5 H T 


X 

-4.36 


+ 10.00 


1 

• 




*WITH CHECK AMPLIFIER 
GAIN = -I 


03 


< 

INT | -0.02 Q l 


* 

-0.21 


-10.00 


: ! 


ii 


04 


1 | 

! _ _ 






i ; 

! j 


05 

06 

_ 


— 




j r - 










| 


07 








v 


i 


08 


i — 

i 






I 


09 




r 






t •; 
i\i 

i S i 


10 


i 

_4 ; _ 






I ! 
j ! 
j | 

; •. i 


... 


1 1 


I ■ ■ 

H - 






T T ™ ■ 

I ! i 

\ i 1 


, 

12 


HG j -Q_/2.5 X I0 4 

i 




- 10.00 


13 




I . : „ 






i j s 


14 


S 


+ 0.02 Q|_ 

L - - - 




+ 10.00 


I i 

i ! •■ 

! | 


ii 


15 


INT 


TIME BASE 

L . . _ _ • ^ 


* 

+ 0.10 


+ 10.00 


T 




16 
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FIGURE 3.7-5. AMPLIFIER ASSIGNMENT SHEET; SURGE TANK 
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TR-48 POTENTIOMETER ASSIGNMENT SHEET 



DATE PROBLEM 



POT 
NO. 


PARAMETER 
DESCRIPTION 


SETTING 
STATIC 
CHECK 


SETTING 
RUN NO.I 


SETTING 
RUN NO. 


SETTING 
RUN NO. 


NOTES 
/Q= 0.1 


POT 
NO. 


00 


Q A /20 A T /3 








► 




00 


.096 — 








01 


Q A /20A T /3 








► 




01 


.096 — 








02 


Q A /20A T £ 


.096 — ' 










02 










03 


0.4 g A u /L/3 


.244 — 






► 




03 










04 


i 


i 1 = 






04 


05 


! 

I 

i. 








05 


06 


i 


1 . - 

i i 
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06 


07 
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»■ 
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1 1 


— .. — I 


. . . . 


. _ ... i _ ' 






11 i 


12 


j 


, • ! . I . _ i 






12 


13 






! 






13 


14 


i 












14 


15 


1= — - ~ t 
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FIGURE 3.7-6. POTENTIOMETER ASSIGNMENT SHEET; SURGE TANK 



SE CT ION IV 
ADVANCED TECHNIQUES 



1. FUNCTION GENERATION 

It is often necessary to generate functions of the independent variable (time) or 
of a dependent variable (voltage) . The following techniques are frequently use- 
ful. 

Analytic Functions 

If a function is given by an analytic expression, such as e^ 21 , sin kx, etc, a 
differential equation can often be found -whose solution is the desired function. 
Mechanizing the differential equation -will then generate the function. 

(1) Generation of Ae"^ 

Let y = Ae-k* 

then y = -kAe"^ = -ky (EQ. 4.1-1) 

The computer circuit is shown in Figure 4*1-1. The function y = Ae ^ is gener- 
ated by including an inverter in the loop. Since this is an expanding exponential 
function, the amplifiers will eventually overload. The function must be scaled 
carefully so that the overload does not occur before the maximum running time, 
which is selected in advance. 

(2) Generation of A sir&fc or A coscot 

Let y = A sir&ot and z = A coscot 

then y = coz and z = -coy (EQ. 4.1-2) 

The computer circuit shown in Figure 4.1-2 arises from the coupled equations 
above or as a solution of the second-order equation y = -to^y. 

(3) Generation of y = where x is a problem variable. The desired 
equation is similar to Section (1) above except that x is a voltage. Since x 
is a function of time; 

y = ke^x = kxy (EQ. 4.1-3) 

The computer circuit is shown in Figure 4. 1-3. Note that the input is x not x. 
This rate input is typical of such systems. In most cases, x is available from 
the circuit that is generating x. 
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(4) Generation of A sin kx where x is a dependent problem variable 



Let y = A sin kx and z = A cos kx 

y * (k A cos kx) x z = (-K A sin kx) x 

y * kzx z = -kyx (EQ. 4.1-4) 

The computer circuit is shown in Figure 4.1-4.. 

Techniques Involved in Using the DFG 

If an analytic expression cannot be found for the function to be generated, then 
use of a DFG is indicated. Since the DFG provides a straight-line segment approx- 
mation to a curve, a certain error is inevitable, and careful programming is nec- 
essary to minimize this error. Often it is not a good idea to generate the desired 
function directly on a DFG. Instead the function should be approximated by a simple 
analytic expression and the DFG used to generate the error or correction term. The 
follcwing examples illustrate the technique. 

Assume the function f (x) in Figure 4.1-5 is to be generated. Since 0 £ x <^ 10 
and 0 £ f(x) £ 75, scale factors of 1 and 1/10 could be used resulting in [x] and 
[1/10 f(x)J as computer variables. 

A DFG could be set up to generate [1/10 f(x)J directly. However, better accuracy 
may be obtained by observing that f(x) has an average value of approximately 70. 
Considering deviations from this mean value, 

Let f^x) = f(x) -70 (EQ. 4.1-5) 

Then -5 1 f-|(x) £ +5, and f*j may be scaled as [2f«j(x)J, Generate [2f-j(x)J on the 
DFG, and then obtain f(x) by solving Equation 4.1^5 for f(x). 

f(x) « f^x) +70 . (EQ. 4.1-6) 

Equation 4. 1-6 must be scaled, just like any other equation. For computer variables 
[1/10 f(x)] and [2 f-j(x)] the scaled equation iss 

[| 0 f(x)J - 1 Q [2F!W] + 7 (EQ. 4.1-7) 

Figure 4.1-6 gives the circuit diagram. Note that any drift or noise occurring in 
generating f-](x.) on the DFG is attenuated, along with f-|(x) itself by a factor of 
20. 

c,o Curve Follower and BIVAR 

A curve follower is available to enable the X-I Plotter to be used as a function 
generator. The curve to be generated is set up by using a fine wire (#18) taped 
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FIGURE 4.1-1. GENERATING THE FUNCTION Ae 




FIGURE 4.1-2. GENERATING THE FUNCTION A sincut 
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FIGURE 4.1-4. CONTINUOUS RESOLUTION CIRCUIT-, 
EQUATION 4,1-4. 
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FIGURE 4.1-5. SAMPLE FUNCTION CURVE 
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FIGURE 4.1-6. COMPUTER DIAGRAM ; EQUATION 4.1-7 



over the desired curve drawn on the plotting surface or by drawing the curve ini- 
tially with conducting ink. The wire or conducting ink curve is energized with a 
high-frequency signal; the pen is replaced by a sensing unit which tracks the curve. 
As the arm is moved back and forth in response to the variable x, the sensing unit 
remains above the curve so that its displacement is f (x) . The pen potentiometer 
is energized with reference voltage, thus the sensing unit output is proportional 
to f (x) • If a variable y is used to energize the pen potentiometer instead of ref- 
erence, the output is proportional to yf (x) , enabling the curve follower to per- 
form multiplication and function generation in one step. 

A BIVAR (Bi-Variant Function Generator) is also available to generate functions 
of two variables. 



2. TRANSFER FUNCTION SIMULATION 



Many engineering problems dealing with linear time-invariant systems are conven- 
iently described in terms of transfer functions. If the transfer functions of a 
system are available, one may readily simulate the system on the TR-48. Two gen- 
eral methods are available; first, the transfer function may be treated as a dif- 
ferential equation and programmed directly using standard amplifiers and potentiom- 
eters; second, passive element networks, containing resistors and capacitors, may 
be used as input and feedback impedances for amplifiers to produce desired trans- 
fer functions. The first method, in general, uses more amplifiers than the second 
but is more convenient to use, especially with unfactored transfer functions. The 
intermediate derivatives are readily available and the introduction of initial con- 
ditions does not cause complications as it does in the second case. On the other 
hand, the use of RC networks affords a saving in the number of amplifiers used. 
Both methods are illustrated here. (The reader should refer to the literature 
cited in the Bibliography for a more comprehensive discussion of transfer function 
simulation.) 

a,o Transfer Function Simulation Using Standard Amplifiers and Potentiometers 

If a transfer function is "treated as a differential equation, it can be programmed 
directly . As an example, consider the transfer function in Equation 4.2-1* 

, e o(p) 1 

G(p)= i^T = 1^ (EQ. 4.2-1) 
Solving for the highest derivative: 

pE = ^ - -° (EQ. 4.2-2) 



Assuming that the highest derivative exists at the input to an integrator, the 
computer circuit shown in Figure 4.2-1 is obtained. 
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Occasionally a system is represented by an unfactored transfer function, such as: 



' P El (p) p3 + ^,2 + E P+ F 4.2-3} 

The following approach is recommended for programming a transfer function of this 
type resulting in the minimum number of components. First, expand the equation 
into differential equation form, giving: 

p% Q +Dp% 0 +EpE 0 +FE 0 « Ap 2 E 1 -+BpE 1 +CB 1 (EQ. 4. 2-4) 

Next, divide the equation by p* 1 "^, -where n equals the highest power of p. (Here 
the divisor is p 2 .) 

pE 0+ DE 0 + E° + == AE., + El + El (EQ. 4.2-5) 

P p 2 P p ^ 

Solving the equation for the highest derivative and grouping terms according to 
powers of p: 

pE Q = AE.j -DE Q + i (BE 1 -EE 0 ) + (CE., -FE Q ) (EQ. 4.2-6) 

Keeping in mind that the notation 1/p represents an integration, the equation can 
then be programmed to yield the computer circuit shown in Figure 4.2-2. 

J2,. Transfer Function Simulation Using RC Networks 

The transfer function of the operational amplifier is: 
E Q Zf 

where Zf and Z^ n are the short-circuit, transfer impedances of the feedback and 
input networks. (The short-circuit transfer impedance of a network is defined as 
the input voltage-output current ratio when the output terminals of the network 
are short-circuited.) Therefore, the relationship expressed in Equation 4.2-7 
can be used to generate transfer functions of the form 

r(r.) - M(p) 



where M(p) and N(p) are the short-circuit transfer impedances of Zf and Z^ n respec- 
tively. A table of short-circuit admittances of some standard networks is in- 
cluded in Appendix 7 and is used in the illustrative example which follows. 
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FIGURE 4.2-2 COMPUTER DIAGRAM •, EQUATION 4.2-6 



A4737 




IUF 



FIGURE 4.2-3 SIMULATION OF TRANSFER FUNCTIONjG (p) = — 

2ll + 0.05p) 



Example; 

Simulate the transfer function 

G(p) - ^°' 1 P) , 
2(1+0.05p) 

The short-circuit admittance of net-work four is given as 

z ss A(1+p6T) 
R-i 

where A = R-j + R 2 , T - R 2 C, © = 



R-j + R 2 



Using this net -work as an input net -work to an amplifier that has Rf as a feedback 
network produces the required transfer function. (See Figure 4.2-3.) 

^o z f R f 

% * " 2 in 3 " A (l+p») 



^f _ 1Q0K = 1 

R<l+R 2 100K+100K 2 

T = R 2 C = (100x10 3 ) (ixlO** 6 ) = 0.1 



©T = — L. (R 2 C) - ZLL - 0.05 
R 1 +R 2 *" 2 



3. REPRESENT AT ION OF DISCONTINUITIES 

Frequently, non-linear effects do not require multipliers and function generators 
for adequate representation on the computer. These effects can be described graph- 
ically by discontinuous straight line input-output relationships. Some of the most 
common discontinuities occurring in physical systems are simple limiting, coulomb 
friction, hysteresis, and gear backlash. All can be represented reasonably veil on 
the computer by standard circuits containing diodes and/ or comparators. Dry fric- 
tion and hysteresis are discussed as typical applications; several other useful 
circuits are shown in the Appendix 8. 

Figure 4.3-1 shows both a diode circuit and a comparator circuit for simulating 
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coulomb friction. 



In equation form, coulomb friction can be represented as: 



f « +c — < 0 

dt 

(EQ. 4.3-U 

f - -o £ > 0 

dt 

The diode circuit (Figure 4.3-1 b) to produce this characteristic is a high- 
gain limited amplifier. Since the gain of the amplifier is very high with no lim- 
iting, the output will very rapidly suing in one direction or the other, depending 
upon the input, until one of the diodes conducts. The diode conduction effectively 
introduces a low resistance in the feedback path of the amplifier, reducing its gain 
very sharply to a very small value. However, the flat parts of the static friction 
characteristic will have a slight rise since the gain of the limited amplifier is 
not exactly zero during the diode conduction because of the finite resistance of 
the potentiometer in the feedback circuit. If this rise in the output characteris- 
tic is considered objectionable, a low value of K (possibly 1) should be used. 

The simulation of static friction using comparators may be accomplished as shown 
in Figure 4.3-1e. If V-j is positive, V 0 = -c. If V-j is negative, V 0 = +c. The 
only error present in the system is in the vicinity of V-j = 0; here V 0 will be 0 
over only a small range, however, this does not usually interfere with satisfactory 
operation in the overall circuit. 



Figure 4.3-2 shows a diode circuit which can be used to represent hysteresis, or 
back-lash . 

This circuit will generate a hysteresis characteristic, but with certain high fre- 
quency limitations. If the input V-j increases from zero, the output V G will re- 
main at zero until the difference between V-j and V 0 is a. When this occurs the 
upper diode will conduct and the integrator will integrate at a rate of 

V-j - V 0 -a 

7~ r„ volts per second, 

[R p + KRjC 

where R p is the forward resistance of the diode, 

R a is the total resistance of the attenuator, 

K is the setting of the attenuator, 

C is the value of integrating capacitor used. 

The circuit in this phase behaves exactly like a simple first-order time lag with 
a time constant [R p + JR a ] C. However, since [R p + KR a ] will normally be very low 
(< 10K C) and C may also be made small, the time constant is very small and can 
usually be considered insignificant. When V^ starts decreasing, the upper diode 
wiH cut off and V 0 will remain constant until an error of -a is developed between 
V-j and V 0 . The lower diode will then conduct causing V c to decrease. The initial 
value of V 0 can be set as an initial condition on the integrator. The integrator 
feedback capacitor should be selected so that the resulting time constant is small. 
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Vo 





+c 




+ v, 




-c 



a. SATURATION OR COULOMB FRICTION 



+ 10 




-10 

b. FEEDBACK LIMITER 




c. RELAY LIMITER 

FIGURE 4.3-1 COMPUTER MECHANIZATION OF SATUR AT ION , OR COULOMB FRICTION 
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i 



+ 10 K^y 




FIGURE 4.3-2 COMPUTER MECHANIZATION OF HYSTERESIS OR BACKLASH 



4. PARTIAL DIFFERENTIAL EQUATIONS 



Partial differential equations can be successfully solved on the analog computer. 
However, a mathematical manipulation of the equation is necessary since the analog 
computer -will integrate -with respect to only one variable, i.e., time. By using 
finite-difference approximations for the derivatives -with respect to all but one 
variable, a set of ordinary differential equations may be -written. 

This technique can be illustrated by using it to solve the equation 

= P^J (EQ. 4.4-1) 

<Jx 2 T ^ t 2 

Equation 4-4-1 describes the motion of an elastic string, where fils the mass per 
unit length of the string, T is the string tension, and y(x,t) is the transverse 
displacement of the string at time t and at a distance x from one end. For such 
an equation, the values of y anddy/d t must be specified for all x at t = 0, and 
the value of y at both ends of the string must be v specified for all t. These are 
the initial conditions and boundary values of the problem respectively. 

Assuming that the initial conditions and boundary values are given, the problem 
can be approached by considering the independent variable x at discrete points 
x Q , x-j, ... , x n , and writing y-? (t) for y(xj_,t) . An ordinary differential equa- 
tion can be -written for each y^(t) by replacing the space derivatives (derivatives 
-with respect to x) by finite-difference approximations. The resulting system of 
ordinary differential equations can be solved on an analog computer. 

To obtain finite-difference expressions for space derivatives, consider the ratio 

y(x+Ax,t) - y(x,t) 
Ax 

This is the slope of the secant line in Figure 4.4-1? and represents the average 
rate of increase of y over the interval from x to x + Ax, for a fixed t. If Ax 
is small, this is an approximation tody^x at (x + Ax/2,t), since x + Ax/2 is the 
midpoint of this interval. A similar approximation holds for the interval from 
x - Ax to x, -whose midpoint is x - Ax/2. Hence the following approximations 
are obtained: 

d7 



^x 



= ?(» * A*,t)-y(x,t) 

(a + &,t) A* 



2 



<^x 



~ 7 U,t) - y(, -Ax,t) M) 



I Ax + \ Ax 
(x- -~,t) 



The accuracy of these approximations depends on the size of Ax. For maximum 
accuracy, Ax should be as small as possible. 
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Similarly, h^j/ £x 2 can be approximated by applying a finite-difference to dy/&x. 
Since x is midway between x - Ax/2 and x + Ax/2: 

"3 r (x,t) ■■ c * * ovx 2 

d x2 ~ (EQ. 4.4-4) 

Replacing the first derivatives by their approximate values given in Equation 4-4-2 
and Equation 4-. 4-3, results in: 

$L%L * xx ~ y(x + A,t) - 2y(x,t) + y(x - Ax,t) 

3 v 2 - ^ ^ TT2 ( E Q. 4.4-5) 

x (Ax) 

This is the desired equation. It expresses the second partial derivative that 
appears in Equation 4.. 4-1 in terms of the values of y(x,t) at discrete points, 
spaced Ax apart. 

Next consider the n + 1 points x 0 , x^, x n , equally spaced so that Ax 

where Ji is the length of the string. Equation 4.4-5 then yiolds: 

\2 

d y±M „ y i+ i - 2yi + yj-i 

(Ax? (EQ. 4.4-6) 

and combining this with Equation 4.4-1 , results in: 

r(Ax)^ 

Since y 0 and y n are given as boundary values no equations are needed for them. 
Suppose the string is stretched taut between two fixed points. Then y 0 « y n = 0 
for all t. These boundary values, together with the system of Equation 4*4-7, de- 
termines the y A for all t. 

Take n = 5. Then Ax S! ^/5. The equations become: 

y<i B §7§ ( y 2 - 2y-i) ^.4-8) 

V 2 * ^4 (73 - 2y 2 + yi ) (EQ. 4.4-9) 

y 3 = 25T (y^ - 2y 3 * y,,) (EQ. 4-4-10) 
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FIGURE 4.4-1 STRING DISPLACEMENT AT A FIXED INSTANT OF TIME 
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FIGURE 4.4-2 COMPUTER DIAGRAM (UNSCALED), EQUATIONS 4.4 -8 THROUGH 4.4-11 



7 U K ^\ (~ 2 U * *3> 4.4-11) 

Thus by consider ing the string displacement at discrete points, a system of 4 
ordinary differential equations is obtained which can be easily solved on the 
computer o The appropriate (unsealed) computer diagram is given in Figure 4-4-2. 

The above method for solving partial differential equations on the analog com- 
puter is the one most frequently used. More information about this technique 
and others may be found in the following Bibliography. 

1o Rogers, A»Eo andT.V. Connolly: "Analog Computation in Engineering Design", 
McGraw-Hill Book Company, Inc., New York 1960 

2. Karplus, "Analog Simulation", McGraw-Hill Book Company, Inc., New 

York 1958 



3. - T -\ckson, AoSoS "Analog Computation", McGraw-Hill Book Company, Inc., 
New York 1960 



5. THE METHOD OF STEEPEST DESCENTS 



The analog computer is basically a device for solving differential equations, 
which makes it especially valuable for the study of dynamic systems. However, 
certain essentially static problems, such as the solution of simultaneous alge- 
braic equations (linear or non-linear), transcendental equations, linear and non- 
linear programming, and curve fitting can also be handled on the analog computer 
by a number of techniques „ One such technique, the method of steepest descents, 
is finding increasing application,, Although the method will be illustrated by 
solving a trivial algebraic system, its scope of application is much more gen- 
er al o 

Consider the set of algebraic equations 

a^ x 1 + a 12 ■ b^ (EQ. 4° 5-1) 



where a^- and a^o &r e both positive- Attempting a solution by means of the boot- 
strap method used for differential equations results ins 



x . b 1 . a 12 , 

l a«j «| a-| -j * 



^ a 22 a 22 1 



>2 - &21 n 



(EQ. 4.5-3) 
(EQ. 4.5-4) 
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A computer circuit for this system is shown in Figure 4.5-1 • If a<|2 and a2-j have 
the same sign, there will be an even number of amplifiers in the loop. A condi- 
tion of positive, rather than negative, feedback will then exist, and the circuit 
will be unstable unless the total loop gain is less than one, that is, unless 

K = <1 (EQ. 4.5-5) 

a 11 a 22 

This requirement is necessary for stability in any loop with an even number of 
summers and no integrators. If it is violated, an overload will quickly arise, 
and no solution will be obtained. Loops containing an odd number of summing 
amplifiers appear to be free of this difficulty. However, this is not necessarily 
true. When no integrator is present in a loop, the phase shift unavoidably pre- 
sent in the amplifiers will produce high-frequency oscillations if the loop gain 
is too high, again preventing the circuit from reaching a stable condition. Thus, 
no matter how many amplifiers are present, one must be very careful to keep the 
loop gain small if a solution is to be obtained. 

For these reasons a naive approach to the solution of algebraic equations is not 
guaranteed of success. A better approach is to program a set of differential 
equations whose steady-state solution satisfies the original algebraic equations. 
In the example given, one can write the differential equations: 

a 11 x 1 + a 12 x 2 - •b 1 + = 0 (EQ. 4.5-6) 

dxo 

a 21 * a 22 * 2 - b 2 ♦ ~4 - 0 (EQ. 4.5-7) 



The circuit for these equations may be obtained from Figure 4*5-1 by substituting 
integrators for the high-gain amplifiers. 

A steady state solution, if one exists, may be found by setting dx-i/dt and dx2/dt 
equal to zero, and such a solution will clearly satisfy Equations 4.5-1 and 4.5-2. 
However, if a-j 2 and a2^ have the same sign, the loop will contain an even number 
of amplifiers, and positive feedback still exists. Again Equation 4.5-5 must 
be satisfied for stability, as may be seen by analytically solving Equations 
4-5-6 and 4.5-7. 

The method of steepest descents is also based on the idea of writing a differen- 
tial equation whose steady-state solution satisfies the algebraic equation. 
However, this method guarantees stability, and the correct solution is obtained 
as a steady-state solution with no restrictions on the coefficients. 

Consider a system of n equations in m unknowns: 

f ± (x v x m ) b ± i - 1, 2, n (EQ. 4.5-8) 



55 



A474I 




FIGURE 4.5-1 COMPUTER DIAGRAM-, EQUATIONS 4.5-3 AND 4.5-4 
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FIGURE 4.5-2 COMPUTER DIAGRAM; EQUATIONS 4.5-17 THROUGH 4.5-20 



It is not necessary to assume that the £^ are linear, nor that m = n. By definition, 
e i (x 1 ... x m ) - f ± (x t ... x m ) - b ± (EQ. 4-. 5-9) 



S (x-, ... x m ) = £ e i ( E Q- 4.5-10) 

i=1 

Note that S(x-j ... x m ) = 0 if and only if the values x-j, . .., x m satisfy Equation 
4.. 5-8. We want to make the Xj ! s vary with time in such a way that lim S— * 0 as 
t-> «>• Since S > 0, S must decrease . which means that dS/dt should be negative. 
From Equation 4.5-10: 



dS ji de 4 

S-aJl H-i (EQ. 4.5-11) 



Since the ej; are functions of x-j, 



de i m ^ £4 dx-? 

Combining Equations 4*5-11 and 4.5-12: 

# = 2£ Si £ (EQ. 4.5-13) 

dt i-1 J-1 ^j dt 

Upon interchanging the order of summation, Equation 4.5-13 becomes: 
^ m dx, G de± 

it = 2 £ dt E £i "5x7 ^-^> 

j 3 1 i=1 

The functions (x-j , ... x m ) and ^ej/3x^ are known, and a set of differential 
equations for the derivatives dx^/dt is desired which will guarantee that dS/dt 
is negative. One way to do this is to write: 

dx_. 11 de* 

-7 1 - -K r s, -3-= (EQ. 4.5-15) 

dt 1 oxj 

i=1 
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as the basic set of differential equations, -where K is a positive constant. Sub- 
stituting 4.5-15 in -4*5-14 yields: 



dt K M dt j 



^ (EQ. 4.5-16) 



■which clearly implies that d£/dt is always negative. Since S is positive and de- 
creasing, it will approach a limit, (positive or zero) -which will be a minimum. 

The method may be summarized as follows: S is an error term (the familiar sum 
of the errors squared) and the vector (x-j ... x m ) moves around in m-dimensional 
coordinate space looking for a minimum S. In the language of vector analysis, 
Equation 4.5-15 says that the point (x ) moves in the negative direction 

along the gradient, of S, so that S decreases as rapidly as possible. Hence, the 
name "steepest descents". The factor K determines the speed with which the vector 
(x-j ... x ) moves. 

If a solution to Equation 4.5-8 exists, then the minimum value of S is clearly 
zero and the x^ will converge to a solution. 

If no solution exists, a best approximation is found (best in the least-square 
sense). If several solutions (or several minima) exist, they will be obtained 
by starting with different initial conditions. There are no problems of instabil- 
ity such as appeared with the other methods. 

Applying this method to the original problem, rewrite Equations 4*5-1 and 4*5-2 
in terms of the errors: 

a^ x^ + a-j2 x 2 ~ ^1 = £ 1 4*5-17) 
a 21 x 1 * a 22 x 2 ~ ^2 = € 2 4*5-13) 

The partial derivatives <^ej_/<^Xj appearing in Equation 4*5-15 are simply the 
coefficients aj_j, and Equation 4*5-15 becomes: 

dx-i 

— L - -K [ £1 ai1 + e 2 a 21 ] (EQ. 4*5-19) 

dxp 

3^ = -K [ C1 a 12 + £2 a 22 ] (EQ. 4*5-20) 

Equations 4.5-17 through 4.5-20 can now be programmed on the computer. The re- 
sulting diagram appears in Figure 4*5-2. 



57 



As pointed out above, the range of applications of this method is much -wider 
than the simple example used as an illustration. Since it is essentially an 
optimization procedure, it is equally applicable to systems with more equations 
than unknowns. Such systems usually have no solution, in vihich case a "solution" 
is found which minimizes the error term S. The following references and those in 
Appendix 9 contain more detailed descriptions of the application of this method to 
least squares curve fitting, linear and nonlinear programming, and similar optimi- 
zation problems. 

1. Rogers, A.E. andT.W, Connolly: "Analog Computation in Engineering De- 

sign", McGraw-Hill Book Company, Inc., New York 19&0 

2. Fifer, Ss "Analog Computation", McGraw-Hill Book Company, Inc. New York 

1961 
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APPENDIX I 
COMPUTER SYMBOLS 



HIGH GAIN AMPLIFIER 



DIODE FUNCTION GENERATOR, 
FIXED AND VARIABLE 



INVERTING AMPLIFIER 



FUNCTION SWITCH 



SUMMING AMPLIFIER 



SOLID-STATE DIODE 



INTEGRATING AMPLIFIER 



■AAA. FIXED RESISTOR 



-O— ^SS N S D | T D P ° TENTI0METER ' ~^\— FACTOR 



Hi 



Lo 



HAND- SET POTENTIOMETER, 
UNGROUNDED 



-0 Y INPUT OF XY PLOTTER 




QUARTER- SQUARE MULTIPLIER, 
B MULTIPLICATION AND DIVISION MODE 




QUARTER- SQUARE MULTIPLIER, 
B SQUARE AND SQUARE- ROOT MODE 



X| X INPUT OF XY PLOTTER 



L RELAY COMPARATOR 



4+i 



'3J 



IN, 



AI-I 



APPENDIX 2 

SIMPLE CIRCUITS USING AMPLIFIERS AND POTENTIOMETERS 



CIRCUIT DESCRIPTION 



CIRCUIT 



COMPUTER SYMBOL 



I. GROUNDED POTENTIOMETER 



■===■ k = POTENTIOMETER 
SETTING 



e O ke 



2. UNGROUNDED POTENTIOMETER 



k 

r 



Hi 



Hi 



e 2 +k (e-e 2 ) 



O — e 2 +k(e,-e 2 ) 

Lo 



ILo 



IOOK 



3. INVERTER 



IOOK 




-e 




-e 



IOOK 



4. MULTIPLICATION BY -10 



IOK 




-ioe 



10 




-ioe 



IOOK 



5. MULTIPLICATION BY - k 
for I < k < 10 



(for k <l use circuit I feeding 
circuit 3 ) 



— SETTING = 



-ke 



ke 



10 



6. MULTIPLICATION BY 2 



IOOK 

IOOK 
e -f-AAAr-i /\. 



-2e 




OR 



e — * 




-2e 



-2e 



IOOK 



7 MULTIPLICATION BY ^ 



IOOK 
IOOK 



lOOKf- *-( >-* — 



e,+e 2 
~2~ 



2 I 



e,+e 2 



a m 
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3. X|x| AND Y|Y| 



YIYJ 
10 




io>r / 



zz 



6-AAArO 

Rf 

O-H-O 

MULTIPLIER 
7.099 



PATCHING DIAGRAM 

IF ONLY ONE ABSOLUTE VALUE CIRCUIT IS 
REQUIRED, ELIMINATE AMPLIFIER 03 AN-B — X~ 
rR- 







U 




-X 2 

QSMI 


B 









(~I0< Y< + I0) 



NOTE 

yIyI _ Y 2 

10 ' 10 

,2 



WHEN Y>0 



- yJyL+X W hen Y<0 



+ X 



-X 

QSMI 




(-I0<X< + I0) 



N0TE xlxL x 2 



10 



-■^ WHEN X>0 



-•*!^=4-^ 2 when X<0 



YM 
10 




+ 10 ,N 



xlxl 

10 



4. SQUARE ROOT 




= Vio|u| 

SEE NOTE 
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X 2 DFG CIRCUITS 




e 0 AMPLIFIER 43 e Q AMPLIFIER 42 




a nr-i 
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3. 



2 

+ - +y 
10 



+ X Y i 

10 



© +IN o 




O- 1 Cm 

X 2 DFG 16.275 



PATCHING DIAGRAM 



•ix 




COMPUTER DIAGRAM 



-I0<X<+I0 





+ 10 


/ 


\ 




s $ 


-10 


-10 


+ 10 



Y POSITIVE 
/ Y NEGATIVE 



NOTE 

X AND Y MUST BE SCALED TO 
PREVENT AMPLIFIER OVER- 
LOAD 



4. 



X 2 
-— -Y 

10 



■ 2-Y 



10 



B, + IN 

-to O 

0, _ R, 



o 



17: 




h&Q -in 

02 



■o 



+,N d 



X 2 DFG 16.275 



tx 




COMPUTER DIAGRAM 

-10 <X< +10 

NOTE 

X AND Y MUST BE 

e„. SCALED TO PREVENT 



PATCHING DIAGRAM 



+ 10 
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AMPLIFIER 

e 




\ + IO ' N 




* Y NEGATIVE 


/ -10 


* Y POSITIVE 
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5. (■*j+Y)FORX<0 



X 2 

— +Y 
10 




0 2 



o 

On 



-IN 
+ IN 



On 
O 



o 
o 

On 



CH 

On 



CM 

X 2 DFG 16.275 



-X 



[ho-y)<x<o] 



IN 



NOTE 

X AND Y MUST BE SCALED TO 
PREVENT AMPLIFIER OVERLOAD 



J(-X-Y)<-lo] 



-10 ^ 



o 

+ 10 



Y POSITIVE 

-- Y NEGATIVE 
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APPENDIX 5 
LOG X AND 4- LOG X DFG CIRCUITS* 



Q. LOG X DFG CIRCUITS 



I. -LOG |0 X FOR +X INPUT AND + LOG |Q X FOR -X INPUT 



•5 LOG ClO | X | ) 



+ X< 



PATCH SAME AS 
AMPLIFIER 35 FOR 
SECOND +LOG X 
CIRCUIT 

PATCH SAME AS 
AMPLIFIER 34 FOR 
SECOND -LOG X 
CIRCUIT 

+ 5 LOG.Ool x|) 



-X 




X-1 IN 

(0<X< + IO) 




-X 



-IN 



(-10^X010) 



-5L0G (10 |X|) 



+ 10 




•10 


\^+IO 


-10 





'IN 




+ 5 LOG (10 |X|) 





+ 10 




+ 10 




-10 



IN 



2. ANTILOG , BASE 10 



-5 LOG (I0|x|) 



X OUT • 




-o ft) 

LOG X DFG 
- 16.276 



-5 LOG (10 J X f ) 
(-10 < LOG <+IO) 




^ ANTILOG^ LOG (10 |Xpi =X 
(0<X <+IO) 
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3. LOG BASE e 



-LOG 



+ X 



.291 



■I- LOG, 



-X 



.35 




+IN 



O 



-10 



-B 

- ,n o ^3 



<3 

+ IN 

o 



R 2 



+10 



-10 



-IN 



O - CM 

LOG X DFG 
16.276 



+ IN 



(0< X < + IO) 




+ 10 


^^flO 


10 




-10 





-IN 



(HO<X<0) 





0 




+ 10 




+ 10 




-10 



921 



2L0G 6 (I0|X|) 



/log BASE 10 

OUT ; SEE 
\EXAMPLE I 




+ 2 L0G o (10 |X|) 

6 



LOG BASE 10 
SEE EXAMPLE I 



4. ANTILOG , BASE e 



X OUT 



.868 

rO"- L0G e 1N 




USE -LOG DFG 
FOR +2 LOGe - 
(10 IXI) INPUT 



.868 

2LOG e (IO |X|)-Q- 




+ ANTILOG [2L0G e (l0|x|)] =X 
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5. MULTIPLICATION (TWO VARIABLES) 




+ X, N ' 



-IOREF 



+ Y IN 



+ IN 


+ LOG X 






R| 



-B. 





LOG X 



-IN 




XY 
10 



6. DIVISION (TWO VARIABLES) 
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+ B_ 
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+ 10 REF 
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-LOG X 





NOTES : 

1. 'M > M 

2. X * 0 
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8. EXPONENTIAL 
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(a) The following Table contains examples of amplifier circuits 
for simulating transfer functions. A more complete listing may 
be found in Jackson, A.S.J "Analog Computation", McGraiHfii|X Book 
Company, Inc., fiev Xork, 1960. 
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(b) The following table contains the short-circuit admittance and component 
values for some useful networks for simulating transfer functions. A more 
extensive listing may be found in Jackson, A.S. : "Analog Computation", and 
Flfer, S.i "Analog Computation". (See Bibliography.) 
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we request your conformance with the following procedure: 
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bers of the major equipment with which the units are to be used, such 
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If at all possible, please include the purchase order or the EAI project 
number under which the equipment was originally procured. 



Your cooperation in supplying the required information will facilitate the processing of your requests and aid in 
assuring that the correct items are supplied. 



It is the policy of Electronic Associates, Inc. to supply equipment patterned as closely as possible to the indi- 
vidual requirements of the individual customer. This is accomplished, without incurring the prohibitive costs 
of custom design, by substituting new components, modifying standard components, etc., wherever necessary to 
expedite conformance with requirements. As a result, this instruction manual, which has basically been writ- 
ton to cover standard equipment, may not entirely concur in its content with the equipment supplied. It is felt, 
however, that a technically qualified person will find the manual a fully adequate guide- in understanding, oper- 
ating, and maintaining the equipment actually supplied. 



Electronic Associates, Inc. reserves 
the right to make changes in design, 
or to make additions to or improve- 
ments in its product without imposing 
any obligation upon itself to install 
them on products previously manufac- 
tured. 
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